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ABSTRACT 

We present the infrared (IR) and X-ray properties of a sample of 33 mid-IR luminous quasars 
> 6 X 10^"^ erg s“^) at redshift 2 ; « 1-3, identified through detailed spectral energy 
distribution analyses of distant star-forming galaxies, using the deepest IR data from Spitzer 
and Herschel in the GOODS-Herschel fields. The aim is to constrain the fraction of obscured, 
and Compton-thick (CT, A^h> 1-5 x 10^^ cm“^) quasars at the peak era of nuclear and star- 
formation activities. Despite being very bright in the mid-IR band, «30% of these quasars 
are not detected in the extremely deep 2 Ms and 4 Ms Chandra X-ray data available in these 
fields. X-ray spectral analysis of the detected sources reveals that the majority (ss67%) are 
obscured by column densities A^h> 10^^ cm“^; this fraction reaches sa80% when including 
the X-ray undetected sources (9 out of 33), which are likely to be the most heavily-obscured, 
CT quasars. We constrain the fraction of CT quasars in our sample to be «24-48%, and 
their space density to be $ = (6.7 ± 2.2) x 10“® Mpc“^. From the investigation of the 
quasar host galaxies in terms of star-formation rates (SFRs) and morphological distortions, 
as a sign of galaxy mergers/interactions, we do not find any direct relation between SFRs 
and quasar luminosity or X-ray obscuration. On the other hand, there is tentative evidence 
that the most heavily-obscured quasars have, on average, more disturbed morphologies than 
the unobscured/moderately-obscured quasar hosts, which preferentially live in undisturbed 
systems. However, the fraction of quasars with disturbed morphology amongst the whole 
sample is «40%, suggesting that galaxy mergers are not the main fuelling mechanism of 
quasars at z « 2. 
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1 INTRODUCTION 

Since the discovery that most, if not all, galaxy spheroids in the lo¬ 
cal Universe host a supermassive black hole (SMBH) whose mass 
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scales with that of the host-galaxy bulge (e.g., 

Magorrian et al.| 

|1998[|Ferrarese & Merritt|2000[[Marconi & Hunt 

2003b it has been 


clear that there is likely to be a tight connection between the evo¬ 
lution of SMBHs and galaxies. Studies of the galaxy population 
across cosmic time have shown that galaxies were more actively 
growing in the past, as the star-formation rate (SFR) density peaks 
at redshift z ~ 2 and rapidly declines to the present day (e.g., 
|Madau et al.|1996l|Hopkins et al.|2006) . A similar evolution is seen 
in the moderate-to-high luminosity active galactic nucleus (AGN) 
population, where most of the SMBHs accrete more rapidly at high 
redshift (z > 1-2), becoming more quiescent at 2 : = 0 (e.g., [Brandt] 
|& Alexander|2015| and references therein). Whether this parallel 
evolution between galaxies and SMBHs is simply due to a larger 
gas supply to feed both black hole accretion and star formation at 
high redshift, or whether there are other processes that self-regulate 
the SMBH and galaxy growth is still uncertain (e.g., [Alexander 
|Hickox|2012[|Kormendy & Ho|201^ . A complete census of AGN 
is crucial to fully understand the accretion history of the SMBHs 
and their role in galaxy evolution. However, as the Universe at 
high redshift is more gas and dust rich, most of the accretion onto 
SMBHs is expected to be heavily obscured, which makes the iden¬ 
tification and study of a large fraction of the AGN population very 
challenging. 

X-ray observations provide arguably the most efficient way 
to trace the growth of SMBHs, as the high energy radiation can 
penetrate large amounts of gas and dust without significant absorp¬ 
tion. Moreover, the contamination from galaxy emission at X-ray 
energies is small as galaxies are typically much fainter than AGN 
in this band. However, when the column density of the obscuring 
material becomes Compton thick (CT; Nn > 1.5 x 10^"^ cm“^) 
even the deepest X-ray surveys often fail to detect the population 
of the most obscured, CT AGN (e.g., |Xue et al.|20T2l see[Br andt &[ 
[Alexander|2015| for a review). A significant fraction of such AGN 
is predicted by all models of the Cosmic X-ray background (CXB) 
to reproduce the observed background spectrum at high energies 
(E « 30 keV; e.g., [Comastri et al.|1995|[Gilli et al.|2007[[Treister[ 
|et al.|200^ , however, the actual contribution from CT AGN is not 
well constrained and varies between the models, due to different 
assumptions, and the space density of CT AGN is still uncertain. 

Since most of the absorbed AGN emission is re-emitted in 
the mid-infrared band (MIR; A « 3-40 fim) by the circumnuclear 
dust heated by the central SMBH, observations in this band can 
potentially identify the elusive CT AGN missed by X-ray surveys 


Stem et al.|2005[|Fiore et al.|2008[ Alexan-] 
der et al.|2008[ |Georgantopoulos et al.|2011] ^ Donley et al.|2008| 


2012[ [Mateos et al.|2013^ . Moreover, since the infrared (IR) band 

is much less affected by extinction than the UV, optical and soft 
X-ray bands, the selection of AGN in the MIR band is almost unbi¬ 
ased against obscured AGN. The downside of using the IR band to 
study AGN is the strong contamination from star formation (SF) in 
the host galaxies, whose cold dust emission tends to dominate the 
MIR to far-IR (FIR, A « 40-250 /rm) bands, peaking at A « 100- 
160 fim. This is especially the case for low-to-moderate luminosity 
AGN, whose IR emission is often swamped by SF. These sources 
are easily missed by typical MIR colour-colour sample selection, 
which are generally used to identify AGN within the star-forming 
galaxy (SFG) popu lation (e.g., [Lacy et al.|2004[ [Stern et al.|2005[ 
[Donley et al. [2012{ [Mateos et al.[2013[. ' ' 

In this paper we investigate a sample of MIR luminous quasars 
(hereafter, “IR quasars”) at redshift 2 « 1-3, i.e., at the peak era 
of AGN and SF activity. We used the deepest available Spitzer 
and Herschel data in the GOODS-Herschel North and South fields 


(hereafter, GH-N and GH-S, respectively; [Elbaz et al.|20lT) to re¬ 
liably identify AGN in the MIR band through detailed spectral en¬ 
ergy distribution (SED) analyses, which allow us to disentangle the 
AGN and SF emissions at MIR and FIR wavelengths. The SED 
analysis is a much more powerful tool than the simple IR colour- 
colour techniques for identifying AGN, as it can recover the hot 
dust emission from the AGN even when it does not strongly dom¬ 
inate the SED (e.g., [Del Moro et al.|[2bl3[ [Rovilos et al.||2014) ; 
however even the SED analyses have limitations and become in¬ 
complete when the AGN component is faint compared to the SE 
component. Using the deepest Chandra X-ray data available across 
the sky (2 Ms in CDF-N; [Alexander et al.[[2003[ and 4 Ms in 
CDF-S; [Xue et al![[2011[ , we investigate the X-ray spectra of the 
IR quasars to constrain the fraction of obscured and CT sources. 
We also investigate the properties of their host galaxies to place 
them in the context of potential SMBH-galaxy evolution scenarios. 
The paper is organised as follows: in Section 2 the data used in 
our analyses are described, while Section 3 details the SED anal¬ 
ysis method and the selection of our IR quasar sample. Section 4 
presents the X-ray spectral analysis and results for our quasars, as 
well as the properties of their host galaxies assessed by the SFRs 
derived from our SED analyses, and their morphologies, obtained 
from visual classification of HST-WFC3 if-band images. In Sec¬ 
tion 5 the results obtained from our analyses are discussed in the 
wider context of the CXB models and the SMBH-galaxy evolu¬ 
tionary scenario. A summary and final remarks are given in Sec¬ 
tion 6. Throughout the paper we assume a cosmological model with 
Ho = 70 km Mpc“^, Dm = 0.27 and Da = 0.73 
[et al.[ |2003[ l. All the errors are quoted at a 90% confidence level, 
unless otherwise stated. 


Spergel 


2 DATA 

2.1 Infrared Observations 


The GOODS fields have been observed in the infrared band by the 
Spitzer and Herschel space observatories, and include the deepest 
IR data currently available. The Spitzer observations jGiavalisco] 
[et al.[ [2004^ cover the MIR band at 3.6, 4.5, 5.8 and 8.0 pm with 
IRAC I Fazio et ^[2004^ and at 24 pm with MIPS ( [Rieke et"^ 
[2004[ > as part of the GOODS Legacy program (PI: M. Dickinson). 
The 24 pm source catalogue was created using a PSF fitting tech¬ 
nique at the 3.6 pm source positions as priors (> 5cr detections). 
For details of the observations and catalogue we refer to|Magnelli[ 
[et al.[p00^ and [Magnelli et al.] ( |2011[ ). In both fields the 24 pm ob¬ 
servations reach a point source sensitivity limit of ~ 30 pJy (5(t). 

FIR observations of the GOODS fields were undertaken by 
Herschel as part of the GOODS-Herschel survey project (GH;[El-[ 


bazet al. 2011 1 . TheGH survey covers the full«10'x 16' GOODS- 


N field (hereafter, GH-N) with PACS ([Poglitsch et al.[2008[l at 100 
and 160 pm and with SPIRE ( [Griffin et al.|2010[ l at 250,'350, 500 
pm, for a total observing time of 124 h. In GOODS-S a smaller re¬ 
gion of « 10 ' X 10 ' at the centre of the field was observed by PACS 
at 100 and 160 pm and the total exposure was 206.3 h (GH-S;[EL[ 
[baz et al.[2014 . We also used the SPIRE 250 pm data observed as 
part of the Herschel Multi-tiered Extragalactic Survey (HerMES; 
[Oliver et al.|2012) across the GH-S field, to better match the data 
available in GH-N. We note that since the SPIRE data suffer from 
significant flux blending issues, we limit our analyses to the 250 pm 
band in both the GH-N and GH-S fields, where these issues can still 
be corrected fairly well. 
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The Sp/fzer-MIPS 24 source positions were used as pri¬ 
ors to estimate the Herschel 100 /rm, 160 and 250 /rm fluxes 
through PSF fitting (courtesy of E. Daddi). In the GH-N field 819 
sources (~42% of the 24 /rm priors) are detected in at least one 
of the Herschel bands down to 5 (j (3a) sensitivity limits of ~1.7 
(~1.2) mjy, ~4.5 (~2.3) mjy and ~6.5 (~4.0) mjy, respectively. 
In the GH-S field there are 591 sources (~34% of the 24 fcm pri¬ 
ors) detected by Herschel (at 100 ^m, 160 fim or 250 /rm); the 
faintest sources are detected down to 5cr (3 (j) flux densities of 1.3 
(0.7) mJy, 2.7 (1.7) mJy, 7.7 (5.0) miy in the three bands, respec¬ 
tively; however, for general sources the flux densities and errors are 
higher, depending on the local source density and the de-blending 
uncertainties. In case of non-detection in any of the Herschel bands, 
appropriate flux density upper limits were calculated for each prior 
position. 

In our analyses we initially considered all the 24 /rm-detected 
sources within the GH fields, namely 1943 sources in GH-N and 
1747 sources in GH-S (for a total of 3690 sources), with or with¬ 
out a Herschel detection. For the sample in GH-N ~9S% of the 
sources have a redshift identification, of which «63% are spec¬ 
troscopic redshifts, while the remaining are photometric (see |Del| 
|Moro et al.|2013| hereafter DM13, for a detailed description of the 
redshift compilation). In the GH-S field, a redshift identification is 
available for ~86% of the sources, of which «95% are spectro¬ 
scopic redshifts compiled from multiple catalogues available in the 


literature ^Le Fevre et al.|2004 

Mignoli et al.|2005[|Vanzella et al. 

2008 

jPopesso et al.||20091 Ba 

estra et al.||2010HSilverman et al. 

2010 

IXia et al.|201 lb and 5% are photometric redshifts. 


2.2 X-ray Observations 


The X-ray counterparts for the 24 ^m detected sources in GH- 
N and GH-S were found by cross-matching the 24 ^m positions 
with the X-ray source catalogues by |Alexander et aL|j2003| l in the 
Chandra Deep Field North (CDF-N), constructed from the 2 Ms 
Chandra observations covering an area of ~448 arcmin^ of the sky, 
and by |Xue et al.| ( |20lT) in the Chandra Deep Field South (CDF- 
S), from the 4 Ms Chandra data over a region of 464.5 arcmin^. 
The CDF-N catalogue contains 503 X-ray detected sources, down 
to flux limits of / 0 . 5-2 kev ~ 2.5 x erg cm“^ 

and / 2 - 8 keV « 1-4 x erg cm“^ The CDF-S X- 

ray source catalogue contains 740 sources down to flux limits 
of /o.5-2keV « 9.1 X 10“^® erg cm“^ and /2-8 keV « 
5.5X 10“^^ erg cm“^ For details on the data reduction and the 
construction of the catalogues we refer to [Alexander et al.|p003| l, 
[Luo et aL] ( |2008| l and jXue et al.| ( |20lT| ). 

For the X-ray detected sources in our sample (see Sect. |3.2| l. 
X-ray spectra have been extracted and analysed. The data were pro¬ 
cessed using the Chandra Interactive Analysis of Observation^^ 
(CIAO; version 4.3 and CALDB 4.4.1.; [Fruscione et al.|2006 tools 
and the ACIS Extract (AE) software packag^ 1 Broos et al. 2010| 
|2012[ >. The Chandra spectra were produced using the AE software 
(version 2011-03-16), extracting source and background spectra in 
each individual observation, as well as the relative response ma¬ 
trices and ancillary files, and combining them using the FTOOLS 
addrmf and addarf. 


^ http://cxc.cfa.harvard.edu/ciao/index.html 

^ The ACIS Extract software package and Users Guide are available at 
http://www.astro.psu.edu/xray/acis/acis_analysis.html 


2.3 Optical/Near-Infrared Observations 

The GOODS-N and GOODS-S fields have been observed by the 
Hubble Space telescope (HST) ACS camera as part of the GOODS 
project ( [Giavalisco et al.j [2004^ , and WFC3 camera as part of the 
Cosmic Near-Infrared Deep Extragalactic Legacy Survey (CAN- 
DELS), a large program consisting of broad-band photometric UV 
to the near-IR imaging of five separate deep extragalactic fields 
(see, [Grogin et al.||20lT] and [Koekemoer et al.||2011| for details 
on the observations and data reduction). The CANDLES HST ob¬ 
servations in the GOODS-N field cover an area of « 10' x 16' 
(CANDELS-Wide) with a deeper smaller region in the centre 
(CANDELS-DEEP; « 6.8' x 10'). The imaging of the GOODS-S 
field covers a similar region as in GOODS-N, but also includes a 
small ultra-deep area in the centre of the field (Hubble Ultra-Deep 
field; 2 ' X 2 . 3 '). In this paper we mainly use data in the F125W 
and F160W bands, which approximately correspond to near-IR J 
and H bands, but also consider data from the optical F606W and 
F850LP bands (corresponding to the V and 2 bands, respectively), 
to corroborate our visual classification of the galaxy morphology 
(Section |4.4[ l. 


3 SPECTRAL ENERGY DISTRIBUTION AND SAMPLE 
SELECTION 

3.1 Infrared SED decomposition 

To constrain the IR emission from the AGN and the host galax¬ 
ies, our SED fitting approach uses an empirical AGN template and 
five different star-forming galaxy templates from [Mullaney et al.j 
( [2011^ , which have been extended down to 3 /rm and to the radio 
band (DM13). The AGN template at A < 40 ^m is represented by 
a broken power law, where the indices have been fixed at the aver¬ 
age values of Fi = 1.7 and r 2 = 0.7 and the p ower-law break at 
A = 19 /rm|^An extinction law ( Draine|2003) is also applied to 
the AGN template, with the extinction parameter Ay free to vary 
in the range of Ay = 0-30 mag jMullaney et al.|20lT] DM13). The 
Spitzer 8, 16 and 24 pm and the Herschel 100, 160 and 250 pm 
photometric points are used to constrain the source SEDs (see Fig¬ 
ures [A^ and |A^ in Appendix); the Spitzer IRAC flux densities at 
3.6, 4.5, and 5.8 pm are not included in the SED analysis because 
these data are likely to be dominated by starlight emission, which 
is not taken into account in our templates. 

The SEDs for all of the 24 ^m-detected sources with redshift 
measurements in the GH-N and GH-S fields have been initially fit¬ 
ted with the SFG templates only, and subsequently with the AGN -F 
SFG templates obtaining a best-fitting solution for each of the SFG 
templates, using a fitting technique (see DM13, for details). We 
adopted the Bayesian information criterion (BIC; |Schwarz|1978) to 
assess the improvement of the fits due to the addition of the AGN 
component to each of the SFG templates, and therefore to identify 


^ We also performed the SED fitting using Fi as a free parameter for the 
sources with enough data points in the MIR band to constrain the power- 
law slope. We note that this is not the case for all the sources, as at 2 > 1.7 
the 8 pm data point falls out of the wavelength range covered by our SED 
templates. Moreover, in some cases the 16 pm data are upper limits and 
therefore cannot provide strong constraints to the MIR power-law slope. 
However, where the fits with Fi as free parameter provided a better repre¬ 
sentation of the data compared to a fixed Fi = 1.7, we adopted these fits 
as the best-fitting solutions. 
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the best-fitting model for each source. Namely, the AGN -|- SFG 
model is assumed as the best-fitting model if for each source: 

(1) ABICi = BICsFG.i — BICAGN-i-SFG,i > 2 (e.g.,|Kass &| 
|Raftery|1995l , where i = 1-5 refers to each of the SFG templates, 
in at least 4 out of the 5 pairs of best-fitting solutions; 

(2) min(xAGN-i-SFG,j) < min(xiFG,i). with j < i, i.e., the 
minimum of the AGN + SFG solutions that satisfy condition 
(1) should be less than the minimum x^ obtained by fitting any of 
the SFG templates alone. 

If any of these criteria is not met we conservatively adopt the simple 
SFG model as the best-fitting model. These criteria are adopted to 
assess the reliability of the AGN component independently from 
the individual SFG templates. In fact, due to the sparse data points 
available for our SED analysis, in some cases, individual fitting 
solutions from different templates yield similarly good fit to the 
data, preventing us from identifying a unique best-fitting solution. 
We refer to Section 3.1 and Appendix A of DM13, for a detailed 
description of the SED fitting approach and tests on its reliability|^ 
After identifying the best-fitting model, we derived the lu¬ 
minosity of the AGN (i.e., corrected for the galaxy contribution) 
at rest-frame 6 fim (Le/jm, agn), the IR luminosity of the AGN 
(T/agn) calculated across the wavelength range A = 8-1000 fim, 
and the IR luminosity of the SFG (A = 8-1000 /rm; hereafter Lsf), 
as weighted means of the values obtained from the five best-fitting 
model solutions (or from all the solutions satisfying condition (1) 
for the AGN -|- SFG model). Where the AGN component was not 
significantly identified an upper limit to the AGN luminosity at 
6 ^m was set as 50% of the total 6 emission (Figure [^, i.e., 
the limit for the SFG component becoming dominant at 6 fim. For 
the cases where the FIR data (100, 160 and 250 /rm) are all up¬ 
per limits and our SED fitting approach did not provide a reliable 
measurement of the SF emission, we estimated an upper limit for 
Lsf by increasing the normalisation of the SF templates until they 
reached any of the FIR upper limits, without exceeding any of the 
MIR data point measurements (including errors). An IR luminos¬ 
ity was calculated from each of the five templates and the highest 
value was then taken as the Lsf upper limit. 

3.2 IR quasar sample 

From the results obtained from our SED decomposition, we se¬ 
lected a sample of luminous AGN with agn > 6 x 

10^^ erg 44.78 in the logarithmic scale; Eig. [^, defined 

hereafter as our “IR quasar” sample. Considering the intrinsic rela¬ 
tion between the 6 fj,m luminosity and the X-ray luminosity ob¬ 
served for unobscured AGN (e.g., |Lutz et al.||2004| |Fiore et al.| 
|200^ [Gandhi et al.|2009| [Mateos et al.|2015| l, our 6 /rm luminos¬ 
ity cut corresponds to an intrinsic X-ray luminosity L 2-10 kev ~ 
2.5 X 10^^ erg s“^, i.e., clearly in the quasar regime even allow¬ 
ing for the ?5i0.4 dex scatter in the relation. We also restricted our 
sample to the redshift range 2 : = 1-3 as we aim to investigate the 
population of luminous quasars at the peak of activity at 2 : ~ 2 
i jMadau et al.|1996| [Hopkins et al.|2006||Ueda et al.|2014[|Brandt] 


^ We note that here we have slightly modified the criteria to evaluate the 
reliability of the AGN component compared to those presented in DM13, 
by adopting the BIG instead of an /-test. The BIG is a more appropriate 
statistical criterion to compare the results from different models than the 
/-test when the models are non-nested. 



Figure 1. Rest-frame 6 fim AGN luminosity, calculated from the SED fits, 
vs. redshift for all of the 24 fim detected sources in the GH fields (~3200 
sources with 2 < 3, out of the total 3690 24 /im-detected sources). The 
black circles represent the sources where we significantly identify an AGN 
component in the IR SED (Sect. |3T}, while the blue open squares represent 
the X-ray detected AGN. For sources where an AGN component was not 
found in the best-fitting SEDs (“IR SFG”), the 6 fim AGN luminosity is 
plotted as an upper limit (grey arrows). The dashed line at log i^Lg^ni = 
44.78 erg s“^ marks our threshold for the IR quasar selection. 

[& Alexander|2015| l. Our sample consists of 33 sources, 19 in the 
GH-N and 14 in the GH-S fields (see Table[T](. The bolometric lu¬ 
minosities (Lboi) of these quasars, estimated as Lboi « 2 x Lagn 
( assuming a bolometric correction of a factor of ~ 8 at 6 ^m, from 
[Richards et al.|2006| and considering the median ratio of ~ 0.3 be¬ 
tween z^Lb/jm, AGN and Lagn estimated for our quasars), span the 
range log Lboi = 45.4—46.5 erg s“^ (see Tablelljl. For all of these 
sources we find that the AGN component largely dominates over 
the SF component at MIR wavelengths, with the AGN contribut¬ 
ing >66% to the total emission at 6 ^m (/agn = 66 — 100%; 
where /agn is defined as the fraction of AGN luminosity over the 
total at rest-frame 6 /im; see Table [TJ. The best-fitting SEDs are 
reported in Appendix [A| (Fig. [A1|[A^ . We note that for 13 out of 
the 33 sources («40%) Spitzer-IRS low-resolution spectroscopy is 
available (A « 3-20 pm; [Pope et al.|2008b| [Murphy et al.|2009| 
[Kirkpatrick et aH2012^ . Although we do not use these spectra to 
constrain the source SEDs, we find very good agreement between 
our SED best-fitting solutions and the Spitzer-IRS spectra, which 
are shown in Fig. [Al[[A2| (Appendix[At, confirming the reliability 
of our SED fitting approach. 

To explore the completeness of the sample we can take into ac¬ 
count that when the AGN is not dominant in the MIR band, i.e., its 
contribution to the total emission is < 50%, our SED fitting pro¬ 
cedure might start failing in identifying the presence of an AGN. 
We therefore looked for all the sources identified as SFGs from 
our SED fitting procedure (i.e., with no significant AGN compo¬ 
nent) with AGN luminosity upper limits above our IR quasar se¬ 
lection threshold of log z^Le^m, agn > 44.78 erg s“^. There is 
only one SFG satisfying this condition, having log z/Le^m, agn < 
44.85 erg s“^. This source is classified as an AGN from its X-ray 
emission (log Lx ~ 43.2 erg s“^). Given its modest X-ray lumi¬ 
nosity and since the Le^jm, agn upper limit is very close to our IR 
quasar luminosity threshold (only a factor ~1.2 higher), the AGN 
within this SFG might still not fulfil our selection criteria, i.e., its 
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Table 1. IR quasar sample 


No 

(1) 

Coordinates 

(2) 

XID 

(3) 

Redshift 

(4) 

-^X 

(5) 

(6) 

/agn 

(7) 

I^AGN 

(8) 

-^bol 

(9) 

(10) 

SFR 

(11) 

Kx.SF 

(12) 

/x,SF 

(13) 

Opt. class 
(14) 

1 

12:36:00.16-1-62:10:47.4 

- 

2.002 s 

<42.30 

45.44 

80% 

46.01 

46.31 

46.21 

723 

41.61 

>20% 


2 

12:36:22.94 4-62:15:26.5 

137N 

2.583 s 

44.80 

45.57 

100% 

46.09 

46.39 

<45.65 

<197 

<41.19 

0% 

BL 

3 

12:36:32.55 4-62:07:59.8 

171N 

1.995 s 

42.69 

45.33 

94% 

45.85 

46.15 

<45.22 

<73 

<40.87 

<2% 

HEG 

4 

12:36:35.60 4-62:14:24.0 

190N 

2.005 s 

42.90 

45.55 

83% 

46.13 

46.43 

46.29 

858 

41.66 

6% 

HEG 

5 

12:36:42.19 4-62:17:11.2 

223N 

2.730 s 

43.93 

44.81 

76% 

45.53 

45.83 

45.70 

221 

41.23 

0% 


6 

12:36:46.72 4-62:14:45.9“ 

243N 

2.004 s 

42.05 

44.90 

74% 

45.44 

45.74 

45.84 

306 

41.33 

19% 


7 

12:36:49.65 4-62:07:38.1 

259N 

1.610 s 

44.04 

45.37 

94% 

45.92 

46.22 

45.51 

145 

41.09 

0% 

HEG 

8 

12:36:55.82 4-62:12:01.2 

287N 

2.737 s 

43.06 

44.96 

69% 

45.66 

45.96 

45.94 

390 

41.41 

2% 


9 

12:36:56.47 4-62:19:37.6 

- 

1.80 p 

<42.33 

45.01 

73% 

45.54 

45.84 

45.95 

396 

41.41 

>19% 


10 

12:36:57.93 4-62:21:28.7 

299N 

2.36 p 

44.12 

44.83 

74% 

45.16 

45.46 

45.88 

332 

41.36 

0% 


11 

12:36:59.32 4-62:18:32.4 

307N 

2.14p 

42.74 

44.96 

83% 

45.68 

45.98 

45.64 

195 

41.18 

3% 


12 

12:37:04.34 4-62:14:46.2 

- 

2.214 s 

<42.09 

45.11 

91% 

45.85 

46.15 

45.41 

114 

41.01 

>4% 

HEG 

13 

12:37:06.87 4-62:17:02.1 

344N 

1.019 s 

44.14 

44.79 

90% 

44.92 

45.22 

45.17 

65 

40.83 

0% 

BL 

14 

12:37:16.67 4-62:17:33.3 

390N 

1.146 s 

43.70 

44.95 

98% 

45.47 

45.77 

44.58 

17 

40.39 

0% 

SE 

15 

12:37:19.87 4-62:09:55.2 

398N 

2.647 s 

43.77 

44.92 

85% 

45.45 

45.75 

45.54 

155 

41.11 

0% 

SE 

16 

12:37:26.50 4-62:20:26.6 

423N 

1.750 s 

42.74 

45.42 

89% 

46.12 

46.42 

45.84 

303 

41.33 

4% 


17 

12:37:39.50 4-62:15:58.6 

- 

2.98 p 

<42.76 

44.84 

66% 

45.57 

45.87 

45.97 

414 

41.43 

>5% 


18 

12:37:42.53 4-62:18:11.8 

459N 

2.309 s 

44.70 

44.98 

93% 

45.25 

45.55 

<45.30 

<88 

<40.93 

0% 

BL 

19 

12:37:57.28 4-62:16:27.4 

478N 

2.922 s 

44.41 

45.01 

66% 

45.48 

45.78 

46.01 

458 

41.46 

0% 

BL 

20 

03:32:09.45 -27:48:06.7 

149S 

2.810 s 

43.96 

45.91 

100% 

46.22 

46.52 

<45.80 

<280 

<41.30 

0% 

BL 

21 

03:32:11.78 -27:46:28.2 

176S 

2.81 p 

43.81 

44.95 

92% 

45.12 

45.42 

<45.39 

<109 

<41.00 

0% 


22 

03:32:20.05 -27:44:47.2 

278S 

1.897 s 

42.43 

44.85 

73% 

45.57 

45.87 

45.82 

291 

41.31 

8% 

HEG 

23 

03:32:23.44 -27:42:55.2 

320S 

2.145 s 

42.24 

44.81 

80% 

45.50 

45.80 

45.59 

174 

41.15 

8% 


24 

03:32:23.72 -27:44:11.8 

- 

2.060 s 

<42.00 

44.82 

83% 

45.38 

45.68 

45.34 

97 

40.96 

>9% 


25 

03:32:24.49 -27:50:45.8 

- 

2.630 s 

<42.06 

44.90 

88% 

45.63 

45.93 

<45.39 

<108 

<40.99 

9% 


26 

03:32:25.14 -27:42:19.2 

344S 

1.617 s 

43.77 

44.89 

66% 

45.44 

45.74 

46.00 

438 

41.44 

0% 

BL 

27 

03:32:25.69 -27:43:05.7 

351S 

2.291 s 

43.79 

45.46 

100% 

45.74 

46.04 

<45.37 

<103 

<40.98 

0% 

HEG 

28 

03:32:28.82 -27:48:29.7 

- 

2.550 s 

<41.95 

45.27 

78% 

45.97 

46.27 

45.93 

375 

41.39 

>28% 


29 

03:32:31.47 -27:46:23.2 

435S 

2.220 s 

42.44 

45.01 

75% 

45.30 

45.60 

45.95 

397 

41.41 

9% 

HEG 

30 

03:32:35.72 -27:49:16.0 

490S 

2.579 s 

42.53 

45.43 

78% 

46.08 

46.38 

46.29 

866 

41.66 

14% 

HEG 

31 

03:32:37.76 -27:52:12.2 

518S 

1.603 s 

44.17 

45.32 

94% 

45.67 

45.97 

<45.61 

<180 

<41.16 

0% 

HEG 

32 

03:32:37.79 -27:42:32.8 

- 

2.975 s 

<42.57 

45.35 

96% 

46.07 

46.37 

<45.78 

<269 

<41.29 

5% 


33 

03:32:49.58 -27:47:14.9 

- 

2.600 s 

<42.33 

44.81 

94% 

45.54 

45.84 

<45.26 

<81 

<40.90 

4% 



NOTES: (1) Source number; (2) MIR 3.6 //m source position, RA and Declination; (3) X-ray source identification number from |Alexander et al.H20Q3) in 
CDF-N and from |Xue et al.|(2011^ in CDF-S; the suffix “N” and “S” refer to North and South, respectively; (4) Redshift; “s” and “p” indicate the 
spectroscopic and photometric redshifts, respectively (Sect. (5) Logarithm of the X-ray luminosity at rest-frame 2-10 keV, extrapolated from the 
measured 0.5-2 keV flux (or upper limit), in units of erg s“^; (6) Logarithm of the rest-frame 6 /im luminosity of the AGN measured from the best-fitting 
SEDs in units of erg s“^; (7) Contribution from the AGN emission to the total (AGN + SF) emission at 6 /im; (8) Logarithm of the AGN IR luminosity 
(8-1000 /im) in units of erg s“^; (9) Logaiithm of the AGN bolometric luminosity, derived from units of erg s“^; (10) Logarithm of the IR 

luminosity (8-1000 /im) of the stai-formation component in units of erg s“^; (11) Star-formation rate in units of Mq yr“^; (12) Logarithm of the rest-frame 
2-10 keV luminosity produced by SF, estimated assuming the|Lehmer et al.||2010| relation between SFR and hard X-ray luminosity; (13) Estimated 
contribution from SF to the rest-frame 2-10 keV luminosity (uncorrected for absorption); (14) Optical spectroscopic classification from |Steidel et al.|^2002 ^, 
[Barger et al.H2003f and |Trouille et al.|f2008^ for the sources in GH-N and from |Szokoly et al.|j2004^ and |Silverman et al^j2Q10| in GH-S: broad-line AGN 
(BL), high-excitation emission line galaxies (HEG); star-forming galaxies (SF). “ For this source we report the X-ray coordinates as the 3.6 /im coordinates 
are centred on a low-redshift galaxy {z = 0.556) lying ps 1.8'' from the X-ray source. Careful inspection of the multi-wavelength images has proved that 
although the Spitz^r-IRAC bands are dominated by the low- 2 : galaxy emission, the emission at A > 8 /im is centred at the X-ray source position rather than 
the 3.6 /im position of the low- 2 : galaxy. We therefore expect the contamination from the low-z galaxy to be very small at the wavelengths used in our SED 

analyses of this IR quasar (see Sect. 


intrinsic power could be below our luminosity cut. For this reason 
we do not include this source in our sample, keeping in mind that 
our selection might not be 100% complete in identifying IR quasars 
when the SF luminosity is very high. 

We note, however, that our IR quasar selection through the 


SED analysis is more reliable than the typical MIR colour-colour 
selections used to identify AGN (e.g., |Stem et al.|2005[|Lacy et al.| 
|2007[ |Donley et al.||201^ , as the SED analysis accounts for the 
source redshifts and corrects for any contamination from SF to 
the MIR emission, contrary to any of the colour-colour selection 
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techniques. Indeed, despite our sources being the most luminous 
quasars at MIR wavelengths within the GH fields, the [Stern et al.| 
( [2005^ AGN wedge would miss ~20% of them, as well as the ma¬ 
jority of the less luminous AGN that we identified from the SED 
analysis (see Fig. [TJ, probably due to the contamination from SF 
and/or starlight to the emission in the MIR bands used in the colour- 
colour plot (typically Spifzer-IRAC bands). Also the [Donley et al.j 
{2Q\2 ) AGN selection criteria would miss several of our IR quasars, 
selecting only ~73% of them. On the other hand, the [Lacy et al.[ 
( [2007^ AGN wedge would select ~97% (all but one) of the IR 
quasars in our sample; however, a large fraction of the sources lying 
within the wedge are not identified as AGN according to our SED 
analysis, nor in the X-ray band, as SFGs at low and high redshifts 
can enter the[ Lacy et al.[ ( (2007| > AGN locus, introducing significant 
contamination from non-AGN (e.g., [Donley et al.|2012^ . 

Of our final sample of 33 luminous IR quasars, 24 (~73%) 
are detected in the X-ray band (Sect. [2^. The X-ray identification 
number of these sources (XID) and their rest-frame 2-10 keV lumi¬ 
nosity (Lx) are reported in Table[^ These luminosities are extrapo¬ 
lated from the observed 0.5-2 keV fluxes reported in the |Alexander[ 
[et al.[ ( [2003^ and [Xue et ak] ^2011^ catalogues for the GH-N and 
GH-S sources, respectively, using a photon index T = 1.8 and ap¬ 
propriate k-correction. This simple extrapolation is performed to 
have a first idea of the discrepancy between the observed X-ray 
and the MIR emission of our quasars. We caution, however, that 
this X-ray luminosity calculation does not account for the possible 
contribution from star formation to the observed X-ray flux (see 
Table[T]and section |43t . Accurate estimates of the observed and in¬ 
trinsic 2-10 keV luminosity for the X-ray detected sources, derived 
from spectral analyses, are obtained in Section [4T| and reported in 
Table 1^ we use these luminosities for further analyses throughout 
the paper. 
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Figure 2. Photon index T vs intrinsic hydrogen column density (Ah) for 
the X-ray detected IR quasars in our sample, obtained using a simple ab¬ 
sorbed power-law model (model 1); the sources XID 243N and XID 320S 
are not included in the plot because the T and Ah parameters could not be 
simultaneously constrained from the spectral fit due to the low numbers of 
counts (< 50 cts; see Sect. [4.1j . Sources in GH-N and GH-S are represented 
by circles and squares, respectively. Red filled symbols indicate sources 
with more than 200 net counts in the i? = 0.5-8 keV spectra; open symbols 
are sources with <200 net counts (i.e., with poor counting statistics). The 
grey shaded region represents the intrinsic photon indices typically found 
for AGN (T = 1.8 it 0. 3; e.g.,[ Nandra & Pounds|1994[ [Mainieri et ^ 

[2002[ [Mateos et al.|2003] [Tozzi et al.|2006| . The quasars with low num¬ 

bers of counts have typically flat photon indices (T < 1.5), suggesting the 
presence of additional spectral components, such as soft scattered emission 
and/or reflection (Sect.[4~i~T). 


4 ANALYSES AND RESULTS 

In this Section we present the X-ray spectral analysis performed 
on all of the X-ray detected IR quasars (24 sources) to constrain 
the absorbing column density (Ah) and therefore derive the in¬ 
trinsic X-ray luminosity. We then estimate the obscured quasar 
(Ah> 10^^ cm~^) and CT quasar (Ah> 1.5 x cm~^) frac¬ 
tions from the X-ray spectral fitting results and from a compar¬ 
ison between the observed (i.e., not corrected for absorption) 2- 
10 keV luminosity and the Lejjm, agn (rest-frame), used as a proxy 
of the intrinsic AGN power. We then investigate in detail the host 
galaxy properties through the SFRs estimated from the Lsf de¬ 
rived from the SED fitting analysis, in relation to the AGN intrinsic 
luminosity and X-ray obscuration. We finally explore the impact of 
galaxy mergers/interactions in quasar hosts by visually classifying 
the host galaxy morphology disturbance of our IR quasars using 
high-resolution HST images. 

4.1 X-ray spectroscopy 

For all of the 24 IR quasars with an X-ray detection in the CDF-N 
and CDF-S catalogues ^Alexander et al.|2003[|Xue et al.|20TT| ) we 
performed X-ray spectral analyses to investigate the high energy 
properties of these sources. For the X-ray undetected sources we 
estimated 3a flux upper limits through aperture photometry on the 
Chandra images, assuming a spectral slope of F = 1.4 for the 
count rate to flux conversion. The 2-10 keV luminosity upper limits 
derived from these fluxes are reported in Table [T] 


The spectra were extracted using the ACIS Extract (AE) soft¬ 
ware package jBroos et al.[[201()l[2012[ l and grouped with a mini¬ 
mum of 20 counts per energy bin in order to adopt ^ fitting statis¬ 
tics. The spectra with limited counting statistics (< 200 net counts 
at 0.5-8.0 keV) were binned with a minimum of one count per bin 
and the Cash statistic ([Cash] |1979J was adopted for the spectral 
analysis with XSPEC (v. 12.8.Ig). 

4.1.1 Basic spectral model fitting 

We initially used a simple absorbed power-law model to charac¬ 
terise the source spectra (model 1), in addition to Galactic ab¬ 
sorption, which was fixed to the average values of An.Cai = 
10^° cm“^ in GH-N and An.Oai = 9.0 x 10^® cm“® in GH-S 
( [Dickey & Lockman[[1990| >. We left the photon index F and the 
intrinsic hydrogen column density (Ah) as free parameters. In Fig¬ 
ure we show the F versus Ah obtained from these fits. For two 
quasars (XID 243N and XID 320S; see Table with very low 
numbers of counts (< 50 net counts, 0.5-8.0 keV), we only fit a 
simple power-law model with no intrinsic absorption as no useful 
constraints could be obtained simultaneously on F and Ah from 
model 1. These two sources are not plotted in Fig.|^ The results 
of the X-ray spectral analysis (including errors at 90% confidence 
level) are reported in Table We note that 14 of the 24 analysed 
sources (~58%) have good counting statistics (net counts > 200; 
see Table [^. For these sources we obtained good spectral con¬ 
straints with our simple model, with spectral slopes consistent with 
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Table 2. X-ray spectral results for all of the X-ray detected sources 


No 

XID 

Net Counts 

F mol 

mol 

Fmo 2 

FVh, mo 2 

/scatt (%) 

log 1/2-10 keV 

log ±2-10 keV.intr 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(V) 

( 8 ) 

(9) 

( 10 ) 

2 

137N 

3587 

1 -yT-l-O.OS 
^•''-0.05 

< 0.8 

- 

- 

- 

44.72 

44.72 

3 

171N 

128 

^•^°-0.46 

<15.3 

1.80 

74 'i+ 6 b.a 

'^••^-32.8 

7.4 ± 1.4 

42.94 

43.79 

4 

I90N 

159 

^■^^-0.74 

20 9 + 1 ^-^ 

1 Q 9 -|-U.t)D 

49.7125-3 

7.9 ±6.1 

43.13 

43.82 

5 

223N 

421 

^•°^-0.30 

9 1 + 6.1 
^•1-5.0 

1 91+0-41 

l-^l-0.34 

10 . 8 li"gi 

7.4 ± 11.6 

43.87 

44.15 

6 

243N 

12 

1 . 861 J-J? 

- 

1.80 

< 8.2 

- 

42.06 

42.11 

7 

259N 

2952 

1.86l°;;3 


1 q':j + U-14 
l-^^-0.13 

10 . 6 i;j 

4.2 ± 1.3 

44.22 

44.49 

8 

287N 

56 

1 01+0.41 
^•^^-0.50 

< 2.6 

- 

- 

- 

43.00 

43.00 

10 

299N 

795 

^•°°- 0.12 

< 1.5 

- 

- 

- 

44.05 

44.05 

11 

307N 

89 

^•^^-0.83 

1 8 + 1 
l9.o_^2.8 

171 -[-1 .cJy 

l-'l-1.29 

39.2l|^:t 

7.2 ±6.5 

42.92 

43.54 

13 

344N 

7007 

^•U9_o.04 

< 0.8 

- 

- 

- 

44.03 

44.03 

14 

390N 

2780 

1 48 + 0-10 

1+-0.4 

l-°^-0.14 

^■^- 1.5 

26.5 ±6.5 

43.89 

44.00 

15 

398N 

414 

1 7f;+0-35 

1 - '^-0.32 

14 5+l’-l’ 
i+-J_6.3 

1 qc+ 0.59 
l-^^-0.48 

99 0+28.9 
ZV.V_;L7.6 

13.4 ± 10.9 

43.75 

44.12 

16 

423N 

114 

0 54+0-56 
^+^-0.41 

< 5.6 

1.80 

1 97 t+ 238.7 
i9/.9_io5 1 

10.9 ± 1.9 

42.92 

43.84 

18 

459N 

3301 

1 76+0-08 
'°-0.06 

< 0.4 

- 

- 

- 

44.62 

44.62 

19 

478N 

890 

2 05+0-^o 

< 1.6 

- 

- 

- 

44.35 

44.35 

20 

149S 

1286 

i-°^-0.l8 

20.6lii 

- 

- 

- 

43.91 

44.37 

21 

176S 

677 

1 00+0.27 
^•^■^-0.25 

13.2^5 ? 

^•^11-0.48 

41.8123-9 

15.1 ± 7.9 

43.70 

44.07 

22 

278S 

165 

^-JJ_1.03 

68-+38:! 

1 - '^- 1.21 

89-81^^;^ 

0.8 ±0.5 

42.89 

44.01 

23 

320S 

48 

fv ^o+O. by 
n./9_o.g7 

- 

1.80 

234.811^1° 

1.5 ±0.2 

41.95 

43.65 

26 

344S 

1680 

2 lO+O-lO 

O+ 0.5 

- 

- 

- 

43.69 

43.73 

27 

35IS 

1378 

1 70+0-18 

1 - '^-0.17 

13.7t|i 

1 76+0-24 

1 - '^^- 0.20 

16.411-6 

3.8 ±4.9 

43.82 

44.15 

29 

435S 

59 

1 1 1+0.46 

1-11-0.44 

< 0.8 

1 89+0-67 

l-o^_0.61 

1100.211235^4 

2.8 ±2.9 

42.46 

44.03 

30 

490S 

82 

U-b/_o,76 

<76.0 

1.80 

2591.412037.2 

0.5 ±0.1 

42.55 

44.86 

31 

518S 

4961 

9 17 -l-U.ue 
^•1 '-0.08 

2 1+0.4 

^•1-0.3 

- 

- 

- 

44.17 

44.27 


NOTES: (1) Source number as in Table^; (2) XID as in Table^ (3) Net counts in the Chandra spectra in the observed 0.5-8 keV energy range; (4) Photon 
index, which was left free to vary in our spectral fits (model 1); (5) Intrinsic hydrogen column density in units of 10^^ cm~^; (6) Photon index obtained from 
model 2; r = 1.80 has been fixed where the spectral fits could not reliably constrain the spectral slopes and Nh simultaneously (Sect. [n~7) ; (7) Hydrogen 
column density obtained from model 2, in units of 10^^ cm~^; (8) Fraction of the intrinsic power-law emission that is scattered at soft energies 
{E <2 keV); (9) Logarithm of the rest-frame 2-10 keV luminosity, in units of erg s“^; (10) Logarithm of the intrinsic, i.e. corrected for absorption, 

2-10 keV luminosity (rest frame), in units of erg s“^. 


the intrinsic values typically found for AGN (F = 1.8 ± 0.3; e.g., 
Nandra & Pounds] 1994[|Mainieri et al.|20Q2[ [Mateos et al.|2Q051 

Tozzi et al. ]2006t . The sources with lower counting statistics (net 
counts < 200) have typically flatter spectral indices, ranging be¬ 
tween r « 0.3-1.4, with large errors. This indicates that the spec¬ 
tra of these sources are harder, and likely more obscured than those 
measured for the sources with higher numbers of counts; however 
the poor counting statistics do not allow us to simultaneously place 
tight constraints on the spectral parameters, T and A^h, as they are 
degenerate for low count spectra (e.g., [Mateos et al.||20()8) . The 
flat r in these sources might also indicate that a simple absorbed 
power-law model does not properly characterise their X-ray emis¬ 
sion and additional spectral components might be required (e.g., 
soft scattered component, or reflection). 

For the sources with a flat photon index (T < 1.5) we there¬ 
fore included an extra power-law component in the model to fit 
the emission at soft energies (F? < 1 keV; hereafter model 2). 
There are several processes that can cause emission at soft X-ray 
energies, such as thermal emission from star forming regions, ra¬ 
dio jets, atomic processes in partially ionised gas, partial covering 
by a neutral absorber, or blurred ionised reflection from the inner 
parts of the accretion disc (e.g., [Turner et al.||1997| [Gierlihski &[ 
[Done|2004> ; however, given the limited numbers of counts avail¬ 


able in our spectra it is not possible to fit more complex models. 
We therefore assume a power-law model as an approximation of 
the soft X-ray emission, with the photon index T tied to that of 
the primary power law (e.g., [Brightman et al.|2014^[Buchner et al.[ 
[2014[ [Lanzuisi et al.[[2015^ to limit the number of free parame¬ 
ters in the model[^ Fitting the spectra with model 2 constrained 
the power-law slopes to more typical values for AGN (see Ta¬ 
ble]^, although the uncertainties on the parameters are still large. 
Where reasonable constraints on F could not be obtained, we fixed 
r = 1.8 to obtain better constraints on the column densities. 
The Wh values resulting from the fitting of model 2 are typically 
higher than those obtained from model 1 and in some cases they 
are very high (Wh> 8 x 10 ^^ cm“^) and possibly consistent with 
Compton-thick absorption. In two cases (XID 435S and XID 490S; 
see Sect. 1^ the column densities estimated by model 2 even 
exceed A'h> 10^^ cm“^. Although these N-n values are not re- 


^ We also tested the spectral fits of model 2 with a soft power-law slope 
free to vary, but no useful constraints could be obtained on F for our 
quasars, due to the limited counting statistics and the high redshift of the 
sources, which shifts most of the soft energy emission {E < 1 — 2 keV, rest 
frame) out of the Chandra sensitivity range. 
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Table 3. Spectral fits of the heavily obscured quasars using PLCABS 
(Yaqoob|1997^ and Torus jBrightman & Nandra|201 1) 




PLCABS 

Torus“ 

No 

XID 

ri> 

Ah° 

F*’ 

Ah° 

3 

171N 

1.8 


1.8 

>332.5 

4 

190N 

i-/O_0.72 

45.0tlti 

1 71+0.71 
^-0.71 

39 r+2’^-6 

11 

307N 

1.8 

^o.o_2g 0 

2 39+0-48 
z.oy_Q 39 

> 603.8“* 

16 

423N 

1.8 

42 7+231.1 
^^•'-28.0 

1.8 

137.4t«5;9 

20 

fogs'* 

1 70+0.16 

21 9+4-0 

1 + 
0 0 

“4 hfk 


22 

278S 

1.8 

84.3ti«;® 

1.8 


23 

320S 

1.8 

181-21^7+ 

1.8 

159.8l^3> 

29 

435S 

^■^^-0.62 

685.4t»«5.o 

1 QQ+0.S4 
^-^^-0.74 

9C1 ^+368.6 

30 

490S 

1.8 

1013.4l«62^" 

1.8 

>429.5 


NOTES: “ To fit the Torus model we fixed the opening angle to 6tor = 60° 
and the inclination angle to flmc = 80°, since these sources have typically 
poor counting statistics and these parameters cannot be constrained from 
the spectral fits. ^ Where sensible constraints could not be obtained simul¬ 
taneously for the photon index F and the column density Nn, we fixed 
r = 1.8. ° The column density N^i is expressed in units of 10^^ cm~^. 

For this source we did not include any soft power-law component as it is 
not needed in the fit. 

liable because the simple parameterisation of the spectra with an 
absorbed power law is not suitable to constrain column densities 
of A^h> few X 10^^ cm“^, they still indicate that these sources 
might indeed be heavily obscured, possibly CT. We investigate this 
hypothesis using more suitable models in Section [4.1.2| In model 2 
we also estimated the fraction of intrinsic emission scattered at soft 
energies (/scatt), calculated from the ratio between the normalisa¬ 
tions of the scattered and intrinsic power-law components, which 
is on average only a few percent, typically <5% (e.g.,|Tumer et al.| 
|1997[ |Ueda et n^|2007^ . However, in some cases we find higher 
scattered fractions (see Table [^, which might be due to the con¬ 
tribution from SF emission to the soft energy spectrum {E < \- 
2 keV). 

We also tested model 2 for sources showing significant X-ray 
absorption (A^h> 10 ^^ cm“^) from the fitting of model 1, but with 
steeper spectral slopes (F > 1.5); the results are reported in Table 
l^for all of the sources where the normalisation of the soft scat¬ 
tered power law could be constrained. Although the column densi¬ 
ties obtained from these fits are typically higher than those obtained 
from model 1, the scatter on the parameters (T and Nn) is larger 
and in all cases the results from the two models are consistent. We 
therefore use the best-fitting solutions from model 1 to estimate the 
luminosity for these sources. 

From the best-fitting models (model 1 for the sources with 
Fmoi > 1-5, and model 2 for the remainders) we calculated 
the rest-frame 2-10 keV luminosity of the X-ray detected IR 
quasars (L2-10 kev), as well as the unabsorbed 2-10 keV luminos¬ 
ity (L2-10 kev.intr), i-S-, cotTected for the best estimate of the Wh 
(Tablel^. 


4.1.2 Spectral models for heavily-obscured AGN 

From the spectral analysis using the simple models 1 and 2, for 
nine sources we obtained very large column density values (see Ta¬ 
ble]^ indicating that they are heavily obscured (defined here as: 
Wh> 2 X 10^^ cm“^), or CT sources (Wh> 1.5 x 10^'* cm“^). 


However, simple absorbed power-law models are not suitable to 
reliably constrain such large column densities, as they do not ac¬ 
curately account for the Compton scattering and reflection of the 
absorbed photons through an optically thick material. We therefore 
modelled the intrinsic emission of these sources using more appro¬ 
priate models, such as: i) PLCABS ( |Yaqoob|1997| >, which simulates 
the attenuation of a power-law continuum by dense, cold matter, 
and ii) Torus ( [Brightman & Nandra|201 1| >, which is a Table model 
based on Monte Carlo simulations of X-ray radiative transfer, self- 
consistently taking into account Compton scattering and iron flu¬ 
orescence emission. To represent the soft X-ray emission we used 
a power-law model, as in model 2. To fit the PLCABS model we 
left Wh and F free to vary, with F tied to that of the soft power-law 
component; where both parameters cannot be constrained due to 
low counting statistics, we fixed F = 1.8. The maximum number of 
scatters in the model was set to 12 (see |Yaqoob|1997| l, and all other 
parameters were fixed to their default values. For the Torus model 
we fixed the inclination angle to an edg e-on view, ^jnc = 80°, and 
an opening angl^of 6tor = 60° (see Brightman & Ueda|2012 1 . 
The Wh and F parameters were left free to vary, with F tied to 
that of the soft scattering component, as for the PLCABS model, or 
fixed to F = 1.8. The best-fitting parameters obtained with these 
two models are reported in Table[^ 

In general there is good agreement between the results ob¬ 
tained from the PLCABS and the Torus models. Although the pa¬ 
rameters are not tightly constrained, as the errors are relatively 
large, all the model solutions resulting from our analyses indicate 
that these quasars are, in fact, heavily obscured and in some cases, 
as for XID 320S, XID 435S and XID 490S, even CT. We note 
that three more sources, XIDs 171N, 307N and 423N are con¬ 
sistent with being CT quasars within the uncertainties on the Wh 
parameter (Table [^. 


4.1.3 Notes on some individual sources 

Although most of the source spectra can be fitted with the simple 
models described in the previous Sections, and sensible results are 
obtained for their spectral parameters (see Tables|^and[^, some of 
the sources show some peculiar characteristics. Here we describe 
in more detail the most interesting sources amongst our X-ray de¬ 
tected IR quasar sample: 

XID 17IN: From the spectral fits using model 2 and PLCABS 
this 2 = 1.995 source shows a very high level of absorption 
8 X 10^® cm“^ (F = 1.8), although not at CT levels. However, the 
residuals in the hard band (E >3 keV, observed frame) suggest the 
presence of a reflection component and an emission line at E ~ 
2 keV (observed frame), typical of CT AGN spectra. Indeed, the 
spectral fit using the Torus model yields a column density lower 
limit which is well within the CT regime. |Georgantopoulos et al.| 
( |2009^ have also classified this source as a CT AGN. The strong 
emission line is consistent with the redshifted iron Ka line (E — 
6.4 keV) with an equivalent width of EW = 1.2lj^;® keV; given 
the limited counting statistics of the spectrum the uncertainties on 
the line EW are large, but still consistent with CT absorption. 

XID 243N: The spectrum of this 2 = 2.004 source has only 
12 net counts (0.5 — 8 keV) therefore a more complex model than 


® We note that we also tried the fits with 9tor = 30°, however, due to the 
limited numbers of counts in the spectra it was not possible to determine 
which of the two geometries provided the best fit. 
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a simple power law cannot be fit to the data. The effective pho¬ 
ton index measured from the spectrum is F = l.SfilJ g®. Fitting 
the spectrum with an absorbed power law with F = 1.8 fixed to 
try to constrain the level of absorption results in an upper limit of 
Nii< 8 X 10^^ cm“^, which is consistent with the source being 
unabsorbed. Flowever, the X-ray luminosity of XID 243N is very 
low compared to the 6 fim luminosity measured from its IR SED 
(see Table [T]) and the ty pical Lx — Le^m relations (e.g., |Lutz et al.| 
|2004[[Flor^t al.|2009[|Mateos et al.|2015^ . The observed X-ray lu¬ 
minosity is consistent with a suppression of a factor of ~ 300 with 
respect to the intrinsic luminosity measured in the MIR (see Table 
1^, and still more than a factor of ~120 considering the intrinsic 
scatter of the Lx — Le^m relation (~ 0.35 dex; e.g.,|Mateos et al.| 
|2015| >. This suggests that XID 243N might be a heavily CT quasar 
(see Sect. |4.2| and Fig.[^ and that the X-ray emission seen by Chan¬ 
dra is only the soft scattered component (typically a few percent of 
the intrinsic nuclear emission; e.g., |Turner et al.|1997[|Ueda et al.| 
|2007[[^ung et al.|2007| and references therein), or some emission 
due to SF in the host galaxy (see Table[T]and Section [43| ). Another 
possibility is that this source is intrinsically X-ray weak, as it has 
been found, for instance, for broad absorption line (BAL) quasars 
(see e.g., |Luo et ar]|2014^ . However, XID 243N has never been 
identified as a BAL quasar to date. Although the hypothesis of in¬ 
trinsically weak X-ray emission cannot be ruled out for this quasar, 
we favour here the heavy-obscuration interpretation (see Sect. |4.2| l. 

XID 287N: Similar to XID 243N, this 2 = 2.737 quasar has 
an apparently unobscured spectrum (with low numbers of counts) 
and an X-ray luminosity that is significantly smaller compared to 
the intrinsic AGN luminosity measured in the MIR band. Assum¬ 
ing a standard value for Lx/Le^jm (e.g., Lutz et ar]|2Q04 Fiore 


|et al.|200^ [Mateos et al.|2015) we estimate that the observed X 
ray luminosity is consistent with being suppressed by a factor of 
~ 35 (see Table |^. This suppression is consistent with the effect 
of CT absorption, but it is less extreme than that observed for XID 
243N. Moreover, no other diagnostics have been found to confirm 
the intrinsic AGN luminosity measured in the MIR band, therefore 
the possibility of intrinsic X-ray weakness cannot be ruled out. 

XID 307N: In the lightly binned spectrum (1 count/bin) of 
this quasar, clear residuals around the observed ~2 keV indicate 
the presence of a strong iron emission line, consistent with being 
Fe Ka at the source redshift of 2 = 2.14. We therefore added a 
Gaussian component to the fit, with E — 6.4 keV and line width 
fixed at CT = 50 eV, to measure the strength of the line. Since 
the redshift of this source is photometric, we initially left the red- 
shift parameter in the model free to vary. From the X-ray spectral 
fit we obtained constraints on the redshift of 2 = 2.18lQ'Qg, in 
good agreement with the photometric estimate reported in Table 
[T] Fixing 2 = 2.18 and fitting the spectrum, however, produces 
a very flat photon index of F = and no constraints 

are found on the Ah parameter. Fixing F = 1.8 yields a col¬ 
umn density of Ah= x cm“^ with some resid¬ 

uals in the hard band (E > 2 keV), suggesting a possible con¬ 
tribution from a reflection component in the spectrum. This hy¬ 
pothesis is supported by the high EW measured for the iron line, 
EW = 1.21q g keV (rest frame), which is typical of CT AGN. 
Indeed, fitting the spectrum with the Torus model results in a col¬ 
umn density lower limit of Ah> 6 x 10^^ cm“^ (see Table [^, 
i.e., in the CT regime. We tested this result also fitting the source 
spectrum with a pure reflection dominated model (i.e., PEXRAV 
model; [Magdziarz & Zdziarski|1995) with the reflection parameter 
R = —1 fixed, plus a Gaussian line, implying column density in 
excess of Ah> 10^® cm“^. This model provides a similarly good 


fit to the data as the previous models, yielding a spectral index of 
F = 2.07 I')q' 43. If 307N is really CT, its intrinsic luminosity 
would therefore be L 2 - 1 Q kev.intr ~ 10^® ergs-F 

XID 390N: This 2 = 1.146 quasar has a relatively strong iron 
Ka line (EW — 0.381 q'44 keV, rest frame, ct = 50 eV fixed), 
which could be broad, as leaving ct free to vary in the fit, we ob¬ 
tain CT < 556 eV. Clear residuals suggest the presence of another 
emission line at L « 5.9 keV (rest-frame) with EW « 0.28 keV, 
which could be a relativistically redshifted iron line from the in¬ 
ner parts of the SMBH accretion disc or an outflow (e.g., [Fabian] 
jet al.|2000| i. Although very interesting, more investigation on the 
emission lines of this quasar is beyond the scopes of this paper. 

XID I49S: This 2 = 2.810 quasar has a column density 
Ah« 2 X 10^® cm“® and a photon index F « 1.7 with all the 
models used to fit the spectrum (i.e., model 1, PLCABS and Torus); 
none of the models required the addition of a soft scattered power 
law. Although the X-ray spectrum does not show signs of any 
emission line, clear residuals around E ~ 1.5-2 keV (E ~ 5.7- 
7.6 keV, rest frame) indicate the presence of a strong absorption 
feature forming a P-Cygni type profile, which indicates the pres¬ 
ence of a high-velocity outflow (e.g., [Pounds & Reeves| |2009[ 
[Tombesi et al.||^10| ). Indeed, this source has been identified in 
the optical band as a BAL quasar jSzokoly et al.[|2004) . A BAL 
quasar interpretation would help explaining why the intrinsic X- 
ray luminosity of this quasar, L 2—10 kev,intr ~ 2 x 10'‘‘‘ erg s’F 
which is corrected for the X-ray absorption derived from the spec¬ 
tral fitting, is lower than that expected from the MIR luminosity 
I'Le/jm « 8 X 10^® erg s“®(see Fig.[^, a s BAL quasars ha ve been 
found to be intrinsically X-ray weak (e.g., [Luo et al. 2014) . p] 

XID 435S: Fitting this source spectrum with model 2 (see Sect. 
4.1 1 results in a column density of Nh— I-IIq 7 ^ 10^^ cm“® 
(with F = 1.89lo g)i which is consistent with heavily CT material. 
Using the PLCABS and Torus models, we obtained smaller Ah 
values of Ah« (3 — 7) x 10®'* cm“® (F « 1.9), but still consistent 
with the source being CT (Table [^. However the spectrum does 
not show a clear signature of a strong Fe Ka line. We note that, 
although these models are too complex for the limited counting 
statistics to provide tight constraints on the spectral parameters, the 
intrinsic 2—10 keV X-ray luminosity corrected for the CT column 
densities derived from the spectral fits is in good agreement with 
the expectations from the MIR 6 fj,m AGN luminosity, according 
to the Lx — Lg^ni relation. XID 435S was also classified as a CT 
AGN also by [Brightman & Ueda[j201^ . 

XID 490S: This 2 = 2.579 quasar shows evidence for a 
strong emission line, consistent with the redshifted Fe Ka line at 
L = 6.4 keV, with an equivalent width oi EW > 1 keV (al¬ 
though the limited counting statistics prevent us from obtaining 
better constraints on the EW). The spectral fitting results from 
model 2 constrained the intrinsic column density largely in the CT 
regime (see Table [^. As in the previous case, fitting the spectrum 
with PLCABS or Torus models constrains the column density to 
lower values than that obtained from model 2 (Ah~ (4 — 10) x 
10®* cm“®), but still consistent with heavy CT columns. Similar 


results were obtained by|Feruglio et al.[([2011[l and [Brightman & [ 
|Ueda| ( |2012[ l, who identified this source as CT. As for XID 435S, 


^ The X-ray weakness of BAL quasars is typically estimated with re¬ 
spect to the optical hand (aoxi e.g., [Luo et al.pOM^ . However, since BAL 
quasars have SEDs consistent with those of typical quasars in the optical- 
to-IR hands (e.g., [Gallagher et al.|2007| , we expect to see X-ray weakness 
relative to the MIR emission too. 
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44.0 44.5 45.0 45.5 46.0 46.5 
log i/L^ (6yU.m; AGN) 

Figure 3. X-ray luminosity (2-10 keV rest-frame) vs. AGN rest-frame 6 fim luminosity (in units of erg s“^) calculated from the SEDs; the observed X-ray 
luminosity (not corrected for absorption) is plotted with open symbols, while the intrinsic luminosity, i.e., corrected for the A^h measured from the X-ray 
spectra (reported in Table[^, is plotted with filled symbols. Unobscured and moderately obscured quasars (A^h < 2 x 10^^ cm“^) are plotted in green, while 
the heavily-obscured sources (A^h > 2 x 10^^ cm~^) are plotted in magenta. Circles and squares con'espond to sources in GH-N and GH-S, respectively. 
The shaded region represents the scatter of the intrinsic relation from|Lutz et al.|^2004[ (dashed line); the dotted line represents the|Fiore et al.| 

j2Q09) relation. 


the intrinsic X-ray luminosity estimated by correcting for the high 
A^h obtained from the spectral fit would place this source on the 
Lx — L 6/im relation (see Fig.|^. 

4.2 Obscured AGN fraction 

From the X-ray spectral analysis of the X-ray detected sources we 
found that a large fraction of the IR quasars (~67%, 16 out of 
24 sources) are absorbed in the X-ray band by column densities 
Ah> 10^^ cm“^. Of these sources, 9 can be classified as heav¬ 
ily obscured (reported in Tablejs}, with Nh> 2 x 10^® cm“^, of 
which three (XIDs 320S, 435S and 490S) have CT column den¬ 
sities, in excess of A^h= 1.5 x 10^^ cm~^ in all the models we 
used (see Sect. |4.l| (; three more sources (XIDs 17IN, 307N and 
423N) are consistent with being CT based on their A^h, depending 
on the adopted model (see Tables and |^, and on the detection 
of a strong iron Ka emission line (EW > 1 keV) in the case of 
XID 171N and 307N. We therefore classify these six sources as 
CT quasars. Due to the limited numbers of counts in the spectra, 
however, the strong Fe Kq line expected in CT spectra cannot be 
detected for all of these CT quasars. 

To assess whether our IR quasars are indeed heavily obscured, 
or CT, we compared the X-ray luminosity (2-10 keV, rest frame) 
with the AGN 6 fim luminosity measured from the best-fit SEDs 
(see Sect. [ 33 . Since the MIR band tends to be much less affected 
by obscuration than the X-ray band, the MIR luminosity is of¬ 
ten used as a proxy for the intrinsic AGN luminosity (e.g., |Lutz| 
let al.|2004| [Alexander et al.|2008||Gandhi et al.|20091 DM13). In¬ 
deed, while unobscured sources are expected to lie on the intrinsic 


Lx-Le^m relation (e.g.. 

Lutz et al.|2004| Fiore et al.|2009||Gandhi| 

|et al.|2009[|Mateos et al. 

20151|Stern|2015|l, heavily obscured and 


CT sources, where the X-ray luminosity is strongly suppressed, are 
expected to have weak observed X-ray emission compared to the 
MIR emission, and therefore lie well below the Lx—Te/jm relation. 
In Figure|^we show the rest-frame L 2 _io keV versus Le^jm.AGN 
for our IR quasars; open symbols represent the observed X-ray lu¬ 
minosity, while the filled symbols represent the intrinsic X-ray lu¬ 
minosity L2_io kev,intr, i.s.. Corrected for the A^h derived from 
the X-ray spectral fitting (see Table|^. The unobscured and moder¬ 
ately obscured sources (shown in green) tend to agree with the local 
intrinsic Lx — Le^m relation, although possibly showing a flatter 
slope than that defined by |Lutz et al.| j2004|l, more in agreement 
with |Fiore et al.|j2009| . All of the heavily-obscured IR quasars lie 
well below the correlation, in the region of the plot expected for 
CT AGN. Correcting their observed X-ray luminosity for absorp¬ 
tion brings these sources much closer to the intrinsic Lx — Le^m 
relation, although they still lie systematically below. It is impor¬ 
tant to note that the A^h of heavily-obscured sources, obtained from 
spectral analyses of i? < 10 keV data alone and with the sim¬ 
ple parameterisation used here (i.e., model 1 and model 2), could 
be underestimated, as shown in several X-ray studies of heavily- 
obscured AGN performed using Chandra or XMM-Newton data to¬ 
gether with hard X-ray data (E > 10 keV) from NuSTAR (e.g., 
Del^^oro^^e£^U|201^ |Gandhi et al.||201^ |Lansbury et al.|[2014[ 

Lansbury et al. |2015> . For the three sources securely identified as 

CT quasars (XIDs 320S, 435S and 490S), the intrinsic X-ray lu¬ 
minosity results in very good agreement with that expected from 
the MIR luminosity; this is also true for XID 17IN and 307N 
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Table 4. Indication of CT quasar classification from various diagnostics for our CT quasars candidates. 


No 

(1) 

XID 

(2) 

Class 

(3) 

Ah 

(4) 

EWpeKa 

(5) 

Lx|6/jm/L 2-10 kEV 
(6) 

-^X|lines 

(7) 

Notes 

(8) 

Refs. 

(9) 

1 

- 


- 

- 

>522 (>233) 



1* 

3 

171N 

CT 

T 

/ 

96 (43) 



2 

6 

243N 

CT 

X 

X 

270 (120) 

45.7 

Lya“ 

3 

9 

- 


- 

- 

>289 0129) 



t 

11 

307N 

CT 

T 


42 (19) 



t 

12 

- 

CT 

- 

- 

>195 (>87) 

44.4 

Lya, C IV, He II, C III]*’ 

I 

16 

423N 

CT 

M2, P, T 

X 

122 (54) 



4* 

17 

- 


- 

- 

>47 (>21) 



t 

23 

320S 

CT 

M2, P, T 

X 

282 (126) 



5t 

24 

- 


- 

- 

>256 0114) 




25 

- 


- 

- 

>273 (>122) 




28 

- 


- 

- 

>820 (>366) 



t 

29 

435S 

CT 

M2, P T 

X 

138 (62) 


He II, C III]“ 

6 

30 

490S 

CT 

M2, P T 

/ 

293 (131) 

44.5 

Lya, C IV, He II, C III], Ne V'* 

6,7* 

32 

- 


- 

- 

>233 (>104) 




33 

- 


- 

- 

> 119 053) 



t 


NOTES: (1) Source number as in Table^ (2) XID as in Table^ (3) CT quasar classification: sources that we consider as secure CT quasars, based on 
various diagnostics explored in our analyses are marked as “CT”; all other sources are still CT quasar candidates; (4) X-ray spectral models yielding A^h 
values consistent, within the errors, with CT absorption A^h> 1-5 x 10^^ cm~^: “M2” refers to model 2, “P” to the PLCABS model and “T” to the Torus 
model (see Sects. |4. 1.1 |and|4. 1 .2[; (5) Sources showing an iron Ka line emission with EW > 1 keV; (6) Factor of suppression of the observed X-ray 
luminosity I/ 2 _io keV compared to the intrinsic X-ray luminosity derived from the 6 ^m luminosity (Ax|6/j.m) using the|Lutz et al.|j2004| relation; the 
values in parentheses correspond to the suppression factors estimated using the lower boundary of the intrinsic Lx — L g relation assuming a scatter of 
~ 0.35 dex; (7) Logarithm of the intrinsic X-ray luminosity derived from optical emission lines, in units of erg s~^; we adopted the average emission-line 
flux ratios from |Netzer et al.| l |2QQ6^ and the [O III];^ 5 oo 7 "X-ray flux ratio from |Mulchaey et al.H1994') ; (8) List of emission lines detected in the optical 
spectra used to infer the intrinsic X-ray luminosity in column (7); (9) Reference of previous work where these sources are classified as CT quasars, or CT 


quasar candidates: 1 =|Alexan(ier et al.H2008K 2 = 

Georgantopoulos et al.H2009|, 3 =|Geor^antopoulos et al.|42011J, 4 =|Georgantopoulos et al. 

|i()ll|; 

5 =[BriEhtman et al.[i2014l, 6 =|BriEhtman & Ueda 

{2012\ and 7 =|Feruglio et al.U201 U. ^ Optical emission line information from|Adams et al.| 

20lU; 

[Alexander et al.[^2008r ‘’[Daddi et al.[j2004j; ““[Szokoly et al.U2004j and[Feruglio et al.[|201 1 

. This source is border-line CT inlBrightman et al. 

12014), 

but consistent with CT absorption within the errors, with A^h= (0.5 — 2.9) x 10^ cm~^. 

These sources also have MIR-to-optical colours typical of 


Dust Obscured Galaxies (DOGs; see Fig. |7]and section[5:^. 


if we correct the X-ray luminosity for the A^h derived from the 
Torus model, which gives logL 2 -io kev,intr ~ 45.4 erg s“^ and 
logL 2 _io kev.intr ~ 45.1 erg s“^, respectively, and possibly for 
XID 423N considering the large scatter on the A^h value. This gives 
an extra indication that these sources are indeed CT. 


We also note that all of the X-ray undetected IR quasars (9 
sources), for which we plot the X-ray luminosity upper limits, lie 
in the region of the Lx — Le^m plane expected for CT sources, 
and therefore can be considered as CT candidates. Two of the X- 
ray undetected quasars, #1 and #12, have been already investigated 
by [Alexander et aL] ( |2008| l and were identified as CT on the basis 
of X-ray analyses and Lx — Le^m analyses. For #12 the intrinsic 
X-ray luminosity has also been estimated from the optical emission 
lines (see Alexander et al.|2008) using the average line ratio from 
[Netzer et al.H2006| l and the relation between the [O IIIJasoot and 
X-ray fluxes from|Mulchaey et al. i 1994 1 , yielding log Lx|iines ~ 
44.4 erg s“^ (Tablej^, in agreement with that estimated from the 
MIR emission. This suggests that this source is intrinsically lu¬ 
minous, further confirming the CT nature of this quasar. The ev¬ 
idence found for #12 also support our assumption that the X-ray 
undetected IR quasars are CT candidates. We therefore infer that 


the fraction of X-ray obscured sources {N\i> 10^ 


) in our 


sample is «76% (25/33 quasars, including the X-ray undetected 
quasars), with Ri55% being heavily-obscured quasars (18/33 with 


Ah> 2 X 10^® cm ^) and 15 out of 33 quasars (« 45%) are likely 
to be Compton thick. 


From the Lx — Le/jm comparison it also emerged that two 
sources having unabsorbed X-ray spectra (XID 243N and XID 
287N; see Table have much weaker X-ray luminosity than that 
expected from their MIR luminosity, although for XID 287N the 
Lx /Lg^ni ratio is not as extreme as for XID 243N (see Fig. [^.XID 
243N was classified as a candidate CT AGN also by|Georgantopou-| 


los et al. 


120111 on the basis of its Lx/Le^m ratio, similarly to our 


analysis. The detection of a strong Lya emission line in the optical 
spectrum, which was attributed to AGN emission rather than SF by 
[Adams et al.| ( [201 l[ l, allowed us to estimate the expected intrinsic 
X-ray luminosity (as for source #12) log Lx pines ~ 45.7 erg s”'^, 
calculated assuming the average AGN Lya to [O III]a5007 line ra¬ 


tio from Netzer et al. 


1 20061 and the [O III]a 5007 to X-ray flux ratio 


from [Mulchaey et a . 1 1994[ (. Although this luminosity is higher 
than that inferred from the MIR emission, possibly due to the un¬ 
certainties on the optical emission line ratios and perhaps to some 
contamination from SF emission to the Lya flux, it supports the 
idea of XID 243N being intrinsically X-ray luminous; we therefore 
consider this source as CT (Table [^. In this case, if we assume 
that the observed E < 8 keV spectrum is only due to scattering 
(~l-3%) of the intrinsic heavily obscured emission, we estimate 
the intrinsic X-ray luminosity to be log Lx « 43.6 erg s“^. This 
is in good agreement with the intrinsic Lx — Le^m relation, sup- 
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Figure 4. Left: SFR vs X-ray column density (A^h) for the IR quasars (derived from the best-fitting of model 1 or model 2); for the X-ray undetected sources, 
which are CT quasar candidates, we fixed JVh at a lower limit of log Nii> 24.3 cm“^. The sources have been divided in three redshift bins: 1.0 < 2: < 1.5 
(green circles), 1.5 < z < 2.5 (blue triangles) and 2.5 < 2 < 3.0 (magenta upside-down triangles). The dashed lines repres ent the main sequence (MS) 
SFRs at redshift 2 = 1, 2 = 2 and 2 = 3 (from bottom to top) for a stellar mass of M* = 9 X 10^® Mq jElbaz et al.||2011|, which is the average stellar 
mass found for X-ray selected AGN at 2 Ri 2 by |Mullaney et al.]j2012^ ; the shaded areas correspond to a factor of 2 above and below the MS at each redshift 
(2 = 1: green, 2 = 2: blue and 2 = 3: pink). Right: SFR vs AGN 6 fim luminosity. Symbols are the same as in Figure]^ The grey dashed line indicates the 
limit where the AGN contributes 50% to the total luminosity at 6 /rm; our IR AGN identification approach will be incomplete for SFRs above this dashed line. 
The vertical dotted line at log = 44.78 erg s“^ marks our IR quasar selection. 


porting our classification of this quasar as CT. For XID 287N no 
other diagnostic has been found to establish whether this source 
is intrinsically X-ray luminous (and therefore likely CT) or X-ray 
weak (e.g., |Luo et ar]|2014^ ; we thus do not include this source 
amongst our CT quasar candidates. 

In Table|^we list all the CT quasars and CT quasar candidates 
identified in our analyses and the various indicators that point to 
a CT quasar classification, i.e., primarily the high column density 
(AiH> 1.5 X 10^'* cm“^) and/or strong iron line (EW > 1 keV) 
derived from the X-ray spectral analysis, the suppression of the ob¬ 
served X-ray luminosity (more than factor of ~ 30 for CT sources) 
compared to the intrinsic X-ray luminosity expected from the 6 ^m 
emission (I/x|6Mm), and other indicators found in the literature, 
such as the X-ray weakness compared to the intrinsic quasar lu¬ 
minosity derived from optical emission lines. The total fraction of 
CT quasars (and CT quasar candidates) in our sample is therefore 
~24-48% (8-16 out of 33 quasars) and the total obscured quasar 
fraction is ^79%. 


4.3 Star-formation rates 


Using the results from our SED fitting analyses (Sect. [33 and 
the X-ray spectral analyses, we explore here the connection be¬ 
tween the quasar properties, such as their luminosity and X-ray 
obscuration, and the host galaxy properties, defined in terms of the 
SFR. From the IR luminosity (8 — 1000 fim) of the SFG com¬ 
ponent, derived from the SED fitting (see Sect. |3.1| l, we calcu¬ 
lated the SFR of the galaxies hosting our IR quasars. We assumed 
a |Salpete r| ^1955| initial mass function (IMF) and used the rela¬ 
tion from Kennicutt| ( |1998| l. The range of SFRs estimated for our 
sources is large, SFR« 17 — 866 Mq yr“^ (see Table^. In Figure 
l^we show the SFR versus the X-ray column density Nu (left) and 
the AGN 6 /im luminosity (right) for the IR quasars in our sample. 


separating the heavily obscured from the unobscured/moderately- 
obscured quasars (magenta and green symbols, respectively; see 
Sect. |4.2| (. In Fig. (left) we divide our sample in three redshift 
bins: 2 = 1.0 —1.5, 2 = 1.5 —2.5 and 2 = 2.5 —3.0 (green circles, 
blue triangles and magenta upside-down triangles, respectively), 
and compare their SFRs with those of typical main sequence (MS) 
SFGs at redshift 2 « 1 — 3 (e.g., |Elbaz et al.|2011[|Mullaney e t al.] 
|2012[ >. In the figure we plot the X-ray undetected sources, which 
are CT quasar candidates, with a column density lower limit fixed 
to log Nu> 24.3 cm~^. On average there is good agreements 
with the SFRs of our IR quasars and the MS in each redshift bin, 
however, at 2 = 1.5 — 2.5 there is a significant fraction of sources 
with SFRs more that a factor of 2 above the MS (6/18 sources in 
this redshift bin, i.e., « 33%), which could be classified as star- 
burst galaxies; these sources also tend have higher X-ray absorp¬ 
tion (Nh> 2 X 10^^ cm~^). On the other hand at 2 = 2.5 — 3.0 
only one (heavily-obscured) quasar out of 13 has SFR typical of 
starburst galaxies at similar redshift, while a significant fraction of 
our IR quasars (5/13 sources, Ri38%) have SFRs below the MS. 
This possibly suggests a bias in our sample, as we are likely to 
miss quasars with very high SFRs (see also Fig.|^ right, and Sect. 
|3.2| (. We caution, however, that given the small number of sources 
in each redshift bin, these results are only tentative. 

Overall, we do not find any significant correlation with the 
intrinsic AGN luminosity, Le/jm.AGN, although the range in lu¬ 
minosity for our IR quasars is relatively narrow; this is in agree¬ 
ment with other studies by |Harrison et al.|p012| > and |Stanley et al.| 
l |2015^ . Interestingly, on average we find that the most obscured 
quasars have higher SFRs than the unobscured quasars (see also 
|Chen et al.|2015^ . However, as the number of sources in our sample 
is small, this difference is mainly driven by three heavily-obscured 
sources with very high SFRs («700-900 Mq yr“^, ~16% of the 
heavily-obscured quasars) that are not representative of the average 
IR quasar host population, for which there seems to be no signifi- 
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cant relation between SFR and X-ray obscuration (e.g. |Lutz et aT] 
|2010| [Rosario et al.|2012[|Rovilos et al.|2012t . 

From the measured SFRs we calculate the expected X-ray 
emission at 2-10 keV (rest frame; I/x,sf) produced by star for¬ 
mation using the [Lehmer et al.j j2010[ l relation between SFR and 
rest-frame 2-10 keV luminosity (equation 4 in their Table 4). The 
values are reported in Tahle[T] We also estimate what fraction of the 
observed X-ray luminosity (at 2-10 keV, rest frame) could there¬ 
fore be due to star formation (see Table[T]). For the majority of the 
sources (~82%) the contrihution from star formation is negligible, 
of the order of <10%. However, for «18% of the sources more 
that 10% of the observed X-ray luminosity could be due to star 
formation and for some of them (#1, 6, 9 and 28) this contrihu¬ 
tion can reach up to «20-30%. All of these quasars are candidates 
CT according to our analyses, of which 3 out of 4 are X-ray unde¬ 
tected. The relatively high contrihution from star formation to the 
2-10 keV luminosity (or upper limit) supports the idea that in these 
quasars we do not detect the intrinsic AGN emission at X-ray en¬ 
ergies, which is likely to he very heavily absorbed, as discussed in 
section |4T3l for XID 243N (#6). 


4.4 Host galaxy interactions 

In order to test the major merger evolutionary scenario for quasars, 
which predicts that most of the accretion onto the SMBH happens 
in heavily obscured environments, due to large amounts of gas 
and dust driven to the centre of galaxies during the merger (e.g., 
[Sanders et al.||1988[|Di Matteo et ^|2005[ [Hopkins et al.[[2006{ 


see also Alexander & Hickox|2012 for a review), we investigate 

here the connection between the X-ray obscuration and the host 
galaxy morphology of our IR quasars. In particular, we aim to ex¬ 
plore whether the most obscured IR quasars are typically found in 
interacting/merging systems. 

We used the high resolution HST images to classify the galax¬ 
ies in terms of disturbance or distortion of their morphology. The 
sources were visually classified mainly using the iT-band images, 
however the supplementary information from the V, z, and J bands 
were also used (see Sect. [2.3[l. Following the visual classifica¬ 
tion method described in Section 4 of |Kocevski et al.[p012[ >, the 
galaxies were initially classified as: 1) “mergers”, where strong 
signatures of distortions and/or multiple nuclei in a single coa¬ 
lescing system were visible; 2) “interactions”, i.e., systems with 
two distinct galaxies clearly showing interaction features, like tidal 
streams; 3) “distorted”, individual galaxies with no visible compan¬ 
ion, but showing asymmetric morphology or distortions; 4) “dou¬ 
ble nuclei”, i.e., galaxies with multiple nuclei in a single coalesced 
system; 5) “close pairs”, where two (or more) galaxies lie within 
~1.5"-2.0" (corresponding to ~13-17 kpc aX z Ki 2), although 
no signs of interactions are visible (which could mean the pairs 
are just line-of-sight alignments); 6) “not-disturbed”, single galax¬ 
ies with no distortion/interaction signatures. We adopted the same 
classification method as in [Kocevski et al.| P012| l in order to al¬ 
low direct comparison between our results and those presented for 
their galaxy and AGN samples, which we consider here as control 
samples for our IR quasars (as they span similar 2 ranges). 

Since the number of sources in our sample, and therefore in 
each class of objects, is small, we grouped the disturbance classes 
listed above in two main groups, following [Kocevski et al.[ ( [2012| l: 

a) “Disturbed", which include all galaxies showing any indication 
of distortion or interaction, i.e., classes 1), 2) 3) and 4); this group 
corresponds to the Disturbed II class in [Kocevski et al.[ \2Q\2) \ 

b) “Undisturbed”, i.e., all sources in classes 5) and 6). Although 


i 


Figure 5. HST F160W {H band) image cutouts (6" X 6") for three of 
our IR quasars showing examples of an “undisturbed” galaxy (#14, left), a 
“companion” class (#5, centre) and a “disturbed” galaxy (#24, right). The 
source numbers are those reported in column (1) in Table^ 



class 5), hereafter defined as “Companion” group, is a sub-class 
of the undisturbed group, we also consider it separately, as despite 
the fact that these close pairs do not show visible signs of distor¬ 
tions/interactions they might still be associated. Examples of these 
three groups are shown in Figure]^ To test the robustness of our 
classification method, we compare our disturbance classification 
for the IR quasars in the GH-S field with the visual classification 
performed by the CANDELS team (e.g., [Kartaltepe et al.[[2012|l. 


vartaltepe 

described in Section 2.1.2 of Rosario et al. j2015 1 jFor details on 
the classification method and the reliability of the classification and 
classifiers we refer to [Kartaltepe et al.[ ( [2012| ). According to the 
CANDELS visual classification, sources are defined as “disturbed” 
if at least 2/3 of the classifiers identified them as such, “interact¬ 
ing” if the interaction metric (IM) is 0.2<IM<0.5 and “mergers” 
if IM>0.5 (see [Rosario et al.|2015^ . We find an agreement of 86% 
between our classification and the CANDELS visual classification 
of the same sources in GH-S; we therefore consider our disturbance 
assessment of the full IR quasar sample as reliable. 

In Figure we plot the fraction of IR quasar host galax¬ 
ies belonging to each of the three disturbance groups described 
above (disturbed, undisturbed and companion). We separated the 
unobscured/moderately-obscured sources (green) and the heavily 
obscured sources (magenta; see Sect. 


4.2 1 . The la errors in the 


plot are calculated following [Gehrels 1 1986[ l. The total fraction 
of IR quasars residing in undisturbed galaxies is « 60%, while 
that of the sources in disturbed/interacting systems is « 40%. 
These fractions are in good agreement with those found by [Ko-[ 
[cevski et al.[j2012) for samples of AGN («56% undisturbed and 
«44% disturbed) and non-AGN hosting galaxies («52% undis¬ 
turbed and «48% disturbed) at similar redshifts to our sample. This 
is true also for the fraction of sources showing a close compan¬ 
ion («12%). Separating the sources between heavily obscured and 
unobscured/moderately-obscured quasars, we find suggestive evi¬ 
dence that the large majority of the unobscured quasars reside in 
undisturbed galaxies (~80%, as opposed to «47% for the heavily- 
obscured quasars), with only ~20% having visible signs of interac¬ 
tions; this fraction seems to be higher («53%) amongst the heavily- 
obscured quasars. Consistent results were found for a sample of 
X-ray selected heavily-obscured AGN in the COSMOS field by 


® The CANDELS visual classification is an ongoing program currently in¬ 
cluding morphological and disturbance classifications for all the galaxies 
with H < 24.5 mag in the GOODS-S and UDS fields, and will soon pro¬ 
vide classifications in all the remaining CANDELS fields. In their catalogue 
every source is inspected by at least 5 classifiers who assigned a number to 
each source as a measure of their disturbance and interaction/merger level 
(or “interaction metric”, IM). 
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Figure 6. Fraction of IR quasar hosts with various disturbance classes, clas¬ 
sified using HST images as described in Sect. |4.4| The fractions and Icr 
uncertainties for the unobscured/moderately-obscured quasar hosts (A^h< 
2 X 10^^ cm~^) are plotted as green triangles and for the heavil y-obs cured 
quasars (A^h^ 2 x 10^^ cm“^) as magenta stai's (see Sect. |4.2| . The 
classes are divided into: “undisturbed”, “companion” (which are a sub¬ 
class of the undisturbed galaxies) and “disturbed”. For comparison the frac¬ 
tions for z ^ 2 AGN and non-AGN samples from |Kocevsld et'n^j2012} 
are also shown (filled black triangles and grey squares, respectively). The 
disturbance fractions for our IR quasars are in good agreement with|Ko-| 
jcevski et al.H201^ . The unobscured/moderately-obscured quasars reside 
preferentially in undisturbed systems, while the heavily-obscured quasars 
are equally found in disturbed and undisturbed systems. 


|Lanzuisi et al. ra- The drop in the number of disturbed galaxies 
amongst the unobscured/moderately-obscured quasars might partly 
be due to the fact that the bright unobscured nucleus might impair 
the detection of faint signs of interactions, such as tidal features; 
we note that particular care has been applied to the classification 
of these sources to minimise this possible bias. Although the errors 
on the fractions of disturbed and undisturbed galaxies are relatively 
large, we find that the difference between unobscured/moderately- 
obscured quasars and heavily-obscured quasars is significant at 
90% confidence level (Fisher exact probability test, P — 0.087). 
However, larger samples are necessary to reliably assess these dif¬ 
ferences at higher confidence levels. 


5 DISCUSSION 

We have investigated the X-ray properties of a sample of 33 MIR 
selected quasars with intrinsic AGN luminosity t/Le/jm > 6 x 
10^^ erg s“ ^ at redshift z « 1 — 3. While at MIR wavelengths these 
quasars are very luminous, the majority have low X-ray luminosi¬ 
ties, and a significant fraction of them (~30%) are undetected in 
the deep Chandra data available in the GH fields, suggesting they 
are heavily obscured in the X-ray band. Indeed, the X-ray spec¬ 
tral analysis and the comparison between the X-ray and MIR lu¬ 
minosities indicate that the vast majority of these IR quasars are 
obscured (up to ~ 79%, 26 out of 33 sources), of which a large 
fraction are candidates to be CT (up to ~ 48%, i.e., 16 out of 33 
quasars). The fractions of obscured and potentially CT quasars are 
much higher than those typically found for optically and X-ray se¬ 
lected quasar samples (see Sect. [ID- We have reliably identified 8 
of these sources as CT from the X-ray spectral analyses (A’h and 


strong Fe Ka line), although some ambiguity still remains for three 
of them (XIDs 17 IN, 307N and 423N) due to the large uncertainties 
on the estimated Wh, from the suppression of the X-ray luminosity 
(a factor of « 30 —1000) compared to the 6 /rm luminosity derived 
from the SED analyses and from the optical emission line luminos¬ 
ity. Further diagnostics are needed for the remaining 8 CT quasar 
candidates to confirm their nature. Eight of our CT quasar candi¬ 
dates have also been classified as CT (or CT candidates) in previ¬ 
ous studies (see Table for reference), however, the remainders 
have never been identified to date. Our results therefore indicate 
that there is a large population of obscured and heavily-obscured 
quasars at 2 « 2 and confirm that the selection of sources in the 
MIR band allows us to find even the most obscured, CT quasars, 
that are missing in the optical and X-ray bands. 


5.1 Obscured quasar fraction: comparison with previons 
studies 


Several studies conducted in the optical and X-ray bands have 
shown that the fraction of obscured AGN decreases with lumi- 


nosity (e.g., |Lawrence & Elvis |1982| |Uec 

a et al. 2003 1 Steffen| 

et al.|2003 |Hasinger|2004| Simpson|2005| 

La Franca et al.|2005 

Treister & Urry 

2005 1 

Hasinger||2008 Uec 

a et al.||2014[. For in- 

stance, [Hasinger 

([2008 

I found that the obscured fraction (/cthin). 


defined as the number of obscured, Compton-thin AGN (Cthin; 
Nn= 10^^ - 10^'‘ cm”®) over the total with Wh< cm”®, 
amongst quasars is more than a factor of 2 lower (/cthin <25% 
at Lx > erg s”^) than that for the less-powerful AGN 

(/cthin ~50% at Lx = 10"^® erg s”^). This trend is also seen 
when selection biases are accounted for (due to the fact that unob¬ 
scured bright quasars are easier to detect out to larger comoving 


volumes, e.g., [Della Ceca et al.|2008^ . The obscured AGN fraction 
seems also to evolve with redshift, becoming larger at high redshift 


(e.g., [La Franca et al. 

2005[[Ballantyne et al.|2006[|Treister & Urry| 

|2006| Hasinger|2008 

Ueda et al.|2014k. The trend of decreasing 


/cthin at high luminosity, however, is found at all redshifts (e.g.. 


|Hasinger|2008||Ueda et al.|2014| ). 

To compare the obscured fraction in our sample with those 
found in previous studies we consider here the Compton-thin 
quasars only, i.e., excluding sources with Nn> 10®^ cm”®. 
Amongst our IR quasars /cthin = 59% — 72%, where the range 
is determined by the lower and upper limit on the number of can¬ 
didate CT quasars in our sample. These fractions are higher than 
those found by [Hasinger] j2008| > for X-ray and optical selected 
Lx « 10'*'^ erg s” quasars at 2 = 1 — 3, /cthin ~ 40% — 55%. 
However, we find good agreement with the jUeda et al.j j2014| l 
estimates in the same redshift range, /cthin ~ 60% — 65% at 
Lx ~ 10"^“^ erg s”^, where they apply a correction for detection 
biases to the observed X-ray absorbed AGN fraction. 

Taking into account the CT sources in our sample the obscured 
quasar fraction reaches « 79% (see Sect. |4.2| l, with CT quasars be¬ 
ing up to ~62% of the total obscured quasars. This fraction is con¬ 
sistent with those found for Seyfert galaxies (Lx < 10'*"^ erg s”^) 
in the local universe, where ~50% of the AGN population shows 
CT absorption (e.g., [Risa li ti et al.||199q [Guainazzi et al.|[2005[ 
[Tozzi et al.||2006] ), possibly suggesting that the relative ratio of 
Cthin and CT sources remains constant across all redshifts and lu¬ 
minosities. Our results also support the assumptions of the CXB 
synthesis models by jGilli et al.jpOOTl l and jUeda et al.|p014| l of an 
equal number of Cthin and CT AGN across all luminosities and red¬ 
shifts. However, we note that different population synthesis models 
(e.g.,|Ballantyne et al.|2006[|Treister et al.|2009't use different as- 
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sumptions on the relative numbers of Cthin and CT AGN (see Table 
6 of |Ueda et al.|20141 . [Buchner et al.|P015| > found that obscured 
AGN (A^h> 10^^ crn“^) at log Lx > 43.2 erg s“^ contribute 
~ 75% of the luminosity density at 2 = 1 — 3, with CT AGN and 
obscured Cthin AGN contributing in equal fractions, while unob¬ 
scured AGN only making up « 25% of the luminosity density. We 
found consistent results for quasars (Lx > erg s“^). 

In the context of galaxy evolution, where the most mas¬ 
sive galaxies, hosts of the most luminous quasars, form earlier 
and evolve more rapidly than the less luminous systems (cosmic 
“downsizing”, e.g., [Cowie et al.||1996[ [Franceschini et al.||1999| >, 
a decreasing fraction of obscured sources with luminosity can be 
explained with the most luminous AGN being more efficient than 
their less luminous analogues in consuming their surrounding gas 
and blowing off their dusty cocoon (through AGN feedback; e.g., 
|Hopkins et al.|[^06[ |Menci et al.||2008] > becoming unobscured, 
bright quasars. This means that the probability of seeing obscured 
high-luminosity quasars is smaller than for lower luminosity AGN, 
as the fast accreting SMBHs spend a shorter time in this “obscured 
phase”. However, at high redshifts we might expect to see more of 
these bright obscured quasars as they have not yet expelled or con¬ 
sumed all their surrounding material. Indeed, [Assef et aL]j2015^ 
found that the space density of very luminous obscured quasars at 
high redshift (2 < a < 4) is comparable to that of unobscured 
quasars of similar luminosities. In this scenario it is perhaps not 
surprising to find such a high fraction of obscured quasars (~ 80%, 
including CT sources) in our sample at z « 2. In fact, there is indi¬ 
cation that CT AGN have a similar evolution with redshift to Cthin 
AGN (e.g., |Brightman & Ueda|20T^ but see also [Buchner et al.| 
|2015| >. The selection in the IR band tends to find the most obscured, 
CT AGN more efficiently than in the optical and X-ray bands, and 
previous studies of MIR selected AGN have indeed found high 
fractions of obscured and CT sources, consistent with our results 


(e.g., Maiolino et al.|2007[|Treister et al.|200^ |Fiore et al.|2008] 

IGeorgantopoulos et al.|201 l[|Stem et al.|2014^ . 


5.2 Comparison with the Dust Obscured Galaxy population 

X-ray obscured AGN have often been found to show extreme MIR- 
to-optical colours (F(24 /im)/F(R); where F(24 /rm) is the flux 
density at 24 /im (observed frame) and F(R) is the flux density 
in the R band), due to strong reddening typically affecting the op¬ 
tical bands (e.g., |Fiore et al.|2009[ [Georganto poulo s et al.|2009| >. 
|Fiore et al.| ( |2008| identified a MIR-to-optical colour selection, 
F(24 pm)/F(R) > 1000, as an efficient way to select heavily ob¬ 
scured, CT AGN, finding a large fraction (~80%) of such sources 
in their sample. Objects with these extreme MIR-to-optical colours 
have often been called Dust Obscured Galaxies (DOGs; e.g., 

|et al.|200^|Pope et al.|200^al ). 

To investigate how our IR quasar selection, and our obscured 
quasar fraction compare to the DOG population, we plot the flux 
density ratio F(24 pm)/F(R) versus the flux density at 24 pm for 
all our IR quasars (Figure]^. The R fluxes for the GH-N sources 
were derived from the Subaru Suprime-Cam R-band magnitudes 
(| Capak et al.[ |2004[ >, while for the GH-S sources they are extrap¬ 
olated from the HST-ACS v and i-band magnitudes jGiavalisc^ 


et al. 112004 1 . We find that a relatively small fraction of 39l 




(the errors on the fraction are estimated from the errors in the 
F(24 pm)/F(R) ratio and Poisson uncertainty) of our IR quasars 
show the extreme colours typical of DOGs. The fraction of heav¬ 
ily obscured quasars amongst them (10/13 sources; « 77%) is in 
good agreement with that found by [Fiore et al.|([^08|l. Of these 
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Figure 7. Mid-infrared to optical flux ratio, F(24 vs the 

Spitzer 24 pm flux density in units of pjy for our IR quasar sample. 
The sources in GH-N are shown as circles and the sources in GH-S as 
squares; heavily-obscured quasars are represented in magenta, while the 
moderately-obscured/unobscured quasars are in green (as in Fig. I^. The 
dashed line marks the flux ratio defining Dust Obscured galaxies (DOGs, 
F(24 pm.)/F{K) > 1000; e.g., [Fiore et al.||2008) . The fraction of IR 
quasars in our sample that can be classified as DOGs is 39^22^®- 
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sources, 8 are amongst our CT quasar candidates (see Table |^. 
However, almost half of the quasars that we identified as heavily 
obscured in our sample would be missed by the DOG colour se¬ 
lection (9/19; « 47%), including «50% of our CT quasar candi¬ 
dates (8/16 sources) I^This indicates that the DOG colour selection 
is highly incomplete in identifying CT quasars and, even though 
some of the most heavily obscured quasars are likely to have ex¬ 
treme F(24 pm)/F(R) colours, this is not true for all of them. 
Indeed, as discussed in section |3.2| the colour selections can be 
strongly affected by contamination from the host galaxy emission 
and by the k-correction, as the source redshift is not accounted for 
in these methods, and therefore cannot provide a clean and com¬ 
plete selection of sources. For our sources we do not find any sig¬ 
nificant relation between the F(24 pm)/F(R) ratio and the red¬ 
shift or the SFRs; further investigation on which of these effects 
have the biggest impact on the DOG colour selection, however, is 
beyond the scopes of this paper. 


5.3 Space density of CT quasars 

Considering the full IR quasar sample, we calculate the sky den¬ 
sity of our sources (z^Le^jm > 6 x erg s~^ai z = 1 — 3) 
within the GH fields (North and South; ~ 260 arcmin^ in total) 
to be «460 deg“^. CT quasars constitute up to « 48% of our 
sources, i.e., up to « 62% of the obscured quasars in our sample 
(26 objects with A^h> 10^^ cm~^). Calculating the space den¬ 
sity of CT quasars (and CT candidates) in our sample (~8-16 
sources) in the comoving volume between 2 ; = 1 — 3 we find 


^ We note that the optical and near-IR photometry of the source #6 (XID 
243N) is contaminated by a nearby low-2: galaxy (see also note in Table^, 
and therefore its F(24 pm)/F(R) ratio, which results to be well below the 
DOG selection threshold, might not be reliable. 


© 0000 RAS, MNRAS 000, 000-000 














































16 A. Del Mow et al. 


(j 

Q. 


cn 

c 

<D 

Q 

<U 

u 

o 

Q. 

to 


10 


10 ' 


-3 


u 

CP 

*> 

o 

E 

o 

o 


10 


-5 . 


10 ' 


' I * ' ' ‘ ‘ ‘ I ' ' ‘ 

“ “ GIlli+07 model 
O FIore+09 
^ Alexonder+08 
□ Tozzi+06 

A Brightmon&Uedo2012 
A This work 


' 




’''9/s 




K> rn^-" ~ ~ - ~ 

' ®"9/s 


0 12 3 4 

Redshift 

Figure 8. Space density of CT quasars at redshift 2 = 1 — 3 estimated 
from our IR quasar sample within the GH fields (black circle). The grey 
dashed lines are the space density predictions from the|Gilli et al.|j2007) 


model for AGN with Lx > 10'*^ erg s ^ 


Lx > 10^® erg s ^andix > 


are also shown: the space density of quasars with Lx > 10^"* erg s ^ 


2 from |Brightman & Ueda (20^ (green triangle), |Fiore et ar]j2009) 
(purple circle)r|AlexanderTr^ j2008 (orange diamond) andlTozzieTar] 
j2006) (blue squares). We al so show the space density of AGN with Lx > 
10^^ erg s“^ al z X 1 from |Tozzi et al.|j2006) . 


<& = (6.7 ± 2.2) X 10“® Mpc~^, where the errors are estimated 
from the lower and upper limits on the number of CT quasars in 
our sample. 

In Figure we plot the comoving CT quasar space density 
estimated in this work and compare it with those found in previ¬ 
ous studies. Our results are consistent with the space density of 
Lx = — 10^® erg s”'^ (intrinsic) CT quasars at 2 « 1.2 — 2.2 

found by |Fiore et al.| (|2009|l, who based their estimates on IR se¬ 
lected sources in the COSMOS field. [Alexander et al.| ( |200^ cal¬ 
culated the volume density of CT quasars at 2 « 2 — 2.5 in the 
CDF-N based on 4 sources identified in the X-ray and IR bands; 
their estimate is somewhat higher than ours, but it is in broad agree¬ 
ment given their large uncertainties. Our estimate is in good agree- 


1 Ms Chandra observations analysed in jTozzi et al.| ( |2006) , while 
the third source was mis-classified as a low-luminosity, unobscured 
AGN. On the other hand, all the CT AGN identified by jTozzi et al.j 
( [2006^ have lower MIR AGN luminosity than our IR quasars (al¬ 
though still consistent with the scatter of the Lx — L^f^m relation), 
and therefore are not selected in our sample, as we are only probing 
the most luminous IR quasars. Two out of the three X-ray detected 
CT quasars we identified in the GH-S field were also classified as 
CT by |Brightman & UedaH201^ , using the same Chandra data as 
in this work; the third source was not included in the Brightman &| 
|Ueda|((^12|l sample, but was analysed by [Brightman et aL j2014 


and classified as a CT quasar. Nevertheless, the overlap between 
our CT quasar candidates and those found by [Brightman & Ueda[ 
\2Q12) is small (« 20 — 25%, only two sources), as most of our 
CT candidates are X-ray undetected, while t he CT quasars (with 
intrinsic Lx > erg s“^) identified by Brightman & Ueda 


( [2012^ have MIR AGN luminosities below our IR quasar thresh¬ 
old (according to our SED analyses). This is perhaps just a selec¬ 
tion effect due to the scatter of the Lx — Te/jm relation; however, 
it might also suggest that the quasar selection approaches in the 
MIR and X-ray bands are complementary, and tend to identify dif¬ 
ferent source populations (e.g., [Brightman & Ueda|2012^ , where 
the MIR selected CT quasars are brighter at MIR wavelengths and 
more obscured (e.g., Ah-J (5 — 10) x 10^"^ cm“^, as the majority 
are X-ray undetected) than the X-ray selected CT quasars, possibly 
because of a larger dust content in the circumnuclear regions. If we 
combine the number of CT quasars identified in our work over the 
GH-S field with those found by [Brightman & Ueda[ ( [2012) within 
the same area (~7 sources within the GH-S area, of which 2 are 
in common with ours), we would obtain a total CT quasar space 
density of <I> = (1.9 ± 0.6) x 10“® Mpc“^. This is in agreement 
with the predictions by |Shi et al.[j2013) based on a joint model of 
the X-ray and IR backgrounds. This suggests that the best estimate 
of the true space density of CT quasars is possibly obtained from a 
combination of the CT populations identified in the IR and in the 
X-ray bands. 


5.4 IR quasars and BH-galaxy co-evolution 

According to some evolution scenarios, the most luminous AGN, 
hosted in the most massive systems, experience a different evo¬ 
lutionary path than that of smaller, less luminous systems. Major 
mergers are predicted to play a significant role in driving gas into 
the centre of the most massive galaxies, causing violent episodes 


2007 1 model for intrinsic Lx > 10'‘‘^ erg s ^ at 2 « 2 and also 

Hopkins et al.[|2006J, especially at high redshifts (e.g., Hopkins[ 

with the values found by|Tozzi et al.|([2006[l and|Brightman & Ueda[ 

et al. 2014|l. On the other hand, secular processes, such as disc 


For a more direct comparison with the X-ray studies of|Tozzi| 
[et al.[ ( |200^ and [Brightman & Ueda[ ( [201^ , we limit our sample 
here to the CT quasars (and candidates) identified in the GH-S field 
(8 sources, three of which are detected in the X-rays). The space 
density estimated by [Tozzi et^(] ( [2006 [ l for CT AGN at 2 « 2, 
is based upon 4 CT AGN identified from X-ray spectral analysis 


in the cenhal 95 arcmin^ in CDF-S (1 Ms). Brightman & Ueda 
( 2012t identified ~ 10 CT AGN in the CDF-S (« 460 arcmin'') 
at 2 « 1 — 3 and intrinsic Lx > 10'*^ erg s“^, yielding a similar 
volume density to [Tozzi et al.[j200^ and to ours. However, it is im¬ 
portant to note that there is no overlap between the CT AGN identi¬ 
fied by [Tozzi et al.| ( [2006) and those identified in this work through 
spectral analysis (i.e., the three X-ray detected sources in GH-S). 
Two of our three X-ray detected CT quasars were undetected in the 


instabilities and galaxy bars, are expected to be the main mech¬ 
anisms driving the gas into the centre of smaller, low-luminosity 
systems, where the SMBH and the galaxy experience a slower and 
more steady growth (see also [Alexander & Hickox|2012[ for a re¬ 
view). To place the IR quasars, which are the brightest MIR sources 
detected in the GH fields, into this evolution scenario, we investi¬ 
gated the host galaxy properties in terms of the SFR and the level 
of disturbance in the galaxy morphology to detect signs of possi¬ 
ble recent interactions/mergers. Our quasars have typically SFRs 
consistent with those of MS SFGs at redshift 2 « 1 — 3, with no 
indication of unobscured quasars having suppressed SF (i.e., below 
the star-forming MS) compared to the most obscured ones (Fig.|^. 
We note, however, that amongst the three quasars with the high¬ 
est SFRs within our sample (SFR^^i 800 Mq yr“^), two are CT 
quasars, one of which is X-ray undetected (#1), and the remainder 
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is heavily obscured by A^h~ 5 x 10^® cm“^. This suggests that 
in these sources large amounts of cold gas over the whole galaxy 
might be partly responsible for the obscuration. However, we do 
not find any direct correlation between the SFRs and A^h over the 
whole sample, suggesting that in general the X-ray obscuration is 
probably confined to the nuclear regions and not related to the gas 
on galaxy scales (e.g., |Rosario et al.|2012[|Rovilos et al.|2012| l. 

As the morphology of the majority of our sources does not 
show any strong distortion features (~ 60%), we can infer that 
mergers are not the main fuelling process amongst the IR quasars 
(Fig. E . This is in good agreement with the merger/interaction frac¬ 
tions found in less luminous AGN and non-AGN hosting galaxies 
at similar redshifts by |Kocevski et al.]p012| >, but it is in contrast 
with the evolution scenario proposed for quasars (e.g., |Di Mat-| 
|teo et al.|2005[ [Hopkins et al.|2006^ . However the fraction of dis¬ 
turbed/interacting systems (~ 40%) is on average higher than those 
found for sources at lower redshifts (« 15 — 20% at z < 1; e.g., 
ICistemas et al.|2011^ , in agreement with the observed increase of 
the galaxy merger rate with redshift (e.g., [Conselice et al.|[2003[ 
|Kartaltepe et al.|2007| ). 

Interestingly, the fraction of unobscured/moderately-obscured 
quasars hosted in disturbed systems is significantly lower («20%) 
than those in undisturbed galaxies (~80%; see Fig. |^, while 
the heavily-obscured quasars are equally found in disturbed and 
undisturbed galaxies, suggesting that the merger fraction in the 
most heavily obscured systems might be higher (see also |Ko-| 
Icevski et al.||2015) . The difference between the two quasar pop¬ 
ulations is significant at the 90% confidence level (P = 0.087, 
see Sect. |4.4^ . If unobscured quasars represent a later stage of 
the SMBH-galaxy evolution compared to the heavily obscured 
quasars, the distortion/interaction features due to mergers or galaxy 
interactions might have faded by the time these quasars are ob¬ 
served as unobscured, as the relaxation time of the galaxy is typ¬ 
ically ~ 200 - 400 Myr (e.g., |Lotzet^L|2^. This might ex¬ 
plain the smaller fraction of disturbed systems we found for the 
unobscured/moderately-obscured quasars. On the other hand, the 
most heavily-obscured quasars might still be at a younger stage of 
evolution after the merger, and therefore their hosts still show the 
signatures of the recent interactions. We note, however, that we did 
not find any direct connection between the SFR and X-ray obscu¬ 
ration in our quasars, and thus no evidence for an enhancement of 
SFR in the heavily obscured quasars in disturbed systems, challeng¬ 
ing the predictions from the quasar evolution scenarios. However, 
as the size of our sample is small, we cannot derive any strong con¬ 
clusion from our results. Further analyses of larger quasar samples 
and comparisons with well selected control samples of galaxies are 
needed to accurately test the SMBH-galaxy evolution scenarios. 


6 SUMMARY AND CONCLUSIONS 

Through detailed SED decomposition analyses of the 24 fim de¬ 
tected sources in the GH-N and GH-S fields we identified a sample 
of 33 luminous IR quasars with AGN MIR luminosity tzLe/jm > 
6 X 10^"* erg s“^ at redshift z ~ 1 — 3. Considering typical in¬ 
trinsic Lx — Lg/jm relations, this luminosity corresponds to an in¬ 
trinsic X-ray luminosity of L 2-10 keV > 2 X erg s \ i.e., 
in the quasar regime. We investigated the X-ray properties of these 
IR quasars through X-ray spectral analysis for the X-ray detected 
sources to accurately derive the obscured quasar fraction and con¬ 
strain the CT quasar population at 2 « 2. Through our SED anal¬ 
yses, which provide reliable measurements of the intrinsic AGN 


luminosity and of the host galaxy SER, and through visual classifi¬ 
cation of the galaxy disturbance morphology using high resolution 
HST images, we also investigated the host galaxy properties in re¬ 
lation to the quasar luminosity and the X-ray obscuration to under¬ 
stand the SMBH-galaxy connection in the context of the evolution 
scenarios. Our main results can be summarised as follows: 

- X-ray spechal analyses of the X-ray detected IR quasars 
(~70% of the sample) show that the majority of the sources (16 
out of 24) are obscured by column densities Wh> 10^^ cm~^, 
with more than half of them (9/16) being heavily obscured (A^h> 
2 X 10^^ cm“^). Given that the fraction of obscured AGN decreases 
at high luminosity (Sect. [ 53 ^ the number of obscured quasars in 
our sample is high; however, we find it is consistent with those 
expected from an increase of the obscured fraction with redshift ( 
after accounting for detection biases, e.g., |Ueda et al.|2014| l. 

- Despite being the brightest quasars at MIR wavelengths de¬ 
tected in the GH fields, about 30% of our IR quasars are not de¬ 
tected in the deepest 2 Ms and 4 Ms Chandra X-ray data. Their X- 
ray luminosity upper limits imply a suppression of the X-ray emis¬ 
sion by factors of ~ 30 —1000 compared to the intrinsic MIR AGN 
luminosity. Erom the comparison of the MIR and observed X-ray 
luminosities for all the IR quasars in our sample we constrain the 
fraction of CT quasars to be ~ 24 — 48% (8-16 quasars, six of 
which are robustly identified through X-ray spectral analyses), i.e., 
up to ~62% of all the obscured quasars in the sample (Sect.|4^. 


- The space density of CT quasars at 2 = 1 — 3 estimated from 
our sample is 4? = (6.7 ± 2.2) x 10“® Mpc“^, which is in gen¬ 
eral agreement with previous results for quasars identified in the 


X-ray or MIR bands (Tozzi et al. 

2006 

|Fiore et al.||2009[|Alexan-| 

|der et ahj |2008| |Brightman & 

Jeda 

20121 as well as with the 


predictions from the population synthesis models (e.g., IGilli et al.| 
|2007| l. However, the small overlap between the CT quasars iden¬ 
tified through our analyses and those found in the X-rays within 


the same fields (e.g., [Brightman & Ueda|[2012^ suggests that the 
space density of CT quasars is likely to be higher than that esti¬ 
mated in this work (or those found in previous studies based on X- 
ray detected sources), and that the combination of X-ray and MIR 
selected CT populations possibly provides the best constraints on 
their true space density in the Universe (Sect. [53) . 

- The SFRs of the IR quasar hosting galaxies are broadly con¬ 
sistent with those of main sequence galaxies at 2 « 1 — 3. In 
general we do not find any direct correlation between the SFR 
and the quasar luminosity or with the X-ray obscuration, with 
only tentative indication for the heavily-obscured quasars (A^h> 
2 X 10^® cm“^) to have higher SFRs than the less obscured ones 
(in particular at 2 = 1.5 — 2.5). The limited number of sources in 
our sample, however, does not allow us to derive any strong con¬ 
clusion from this result (Sect. |43) . 

- From visual classification of the galaxy morphology to iden¬ 
tify signs of interactions/mergers we find that ~40% of our quasars 
have a disturbed morphology, while ~60% do not show any sign 
of distortions, indicating that mergers are not the main mech¬ 
anism fuelling the SMBH in luminous quasars. This is consis¬ 
tent with the results found for lower-luminosity AGN (e.g. |Ko-| 
jcevski et al.|2012) . However, the interaction/merger fraction seems 
to be lower for unobscured/moderately-obscured quasars than for 
heavily-obscured quasars, possibly suggesting that large amounts 
of gas and dust in the galaxy, as a result of recent merger or galaxy 
interactions, are partly responsible for the high level of obscura¬ 
tion in these quasars, while the unobscured/moderately-obscured 
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quasars might be at a later evolutionary stage, where the signatures 
of galaxy interactions/mergers have already faded (Sect. [4!4) . 

Further progress can be made in the near future to constrain the 
CT quasar population using the extremely deep 7 Ms Chandra data 
recently obtained in the CDF-S (PI. W. Brandt). These data will 
allow us to confirm the CT nature of the CT quasar candidates we 
identified in the GH-S field (e.g., providing better spectra and de¬ 
tecting some of the X-ray undetected systems), further validating 
our analyses and results. Moreover, the advent of the new X-ray 
and MIR observatories, such as eROSlTA, ATHENA and the 
James Webb Space Telescope (JWST) will provide unprecedentedly 
large samples of high -2 quasars to allow studies of MIR and X-ray 
selected CT quasars, placing tighter constraints of their space den¬ 
sity and on the different properties characterising sources selected 
in the MIR versus those selected in the X-ray bands. 
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APPENDIX A: INFRARED SPECTRAL ENERGY 
DISTRIBUTIONS OF THE IR QUASARS 
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Figure Al. Spectral energy distribution of the IR quasars in GH-N. Filled circles are the Spitzer 8.0, 16 and 24 pm and Herschel 100, 160 and 250 pm used 
to constrain the SEDs. The Spitzer-\KAC 3.6, 4.5 and 5.8 pm data (open triangles) and the Herschel 350 and 500 pm (open circles) are also shown, although 
these data are not used in the SED fitting process (Sect. ED- For 8 out of 19 IR quasars in GH-N (~42%) Spitzer-IRS spectra are available and they are shown 
here (cyan line) to demonstrate the agreement with the best-fitting SEDs. The dotted curves and the dashed curves represent the AGN and SFG templates, 
respectively, while the solid curve is the total best-fitting SED. The purple stars, which are only shown for the X-ray detected IR quasars, represent the 6 pm 
flux we would expect from the observed X-ray luminosity adopting the|Lutz et al.|J2004| relation (for XID 243N the star falls below the plotted flux range). 
On the top left corner of each plot we report the source number and the redshift, while on the top right comer there is the XID from | Alexander et al. 1^2003*) for 
the X-ray detected quasars. 
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Figure A2. Spectral energy distribution of the IR quasars in GH-S. Symbols are the same as in Fig. 1^ The 5p//z^r-IRS spectra, which are available for 4 out 
of the 14 IR quasars in GH-S (~29%) are shown (cyan lines) to demonstrate the agreement with the best-fitting SEDs. For XIDs 320S, 435S and 490S, which 
we robustly identified as CT quasars, the purple stars, which represent the 6 /im flux we would expect from the observed X-ray luminosity adopting the |Lutz| 
|et al.|^2004[ relation, falls below the plotted flux range. On the top left corner of each plot we report the source number and the redshift, while on the top right 
comer there is the XID fromjXue et al.|(2011) for the X-ray detected quasars. 
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